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Abstract 


This  thesis  investigates  enhancements  of  an  existing  numerical  predictor-corrector 
aerocapture  guidance  algorithm  (PredGuid).  The  study  includes  implementation  of  an  energy 
management  phase  prior  to  targeting  with  a  generic  method  of  transition  and  replacement  of 
heuristic  features  with  more  generic  features. 

The  vehicle  response  during  energy  management  was  modeled  as  a  second-order 
spring/mass/damper  system.  Phase  change  occurred  when  two  conditions  were  met:  First,  the 
vehicle  could  fly  a  constant  bank  angle  of  1 10°  for  the  remainder  of  the  trajectory  and  have  the 
resulting  apogee  below  or  within  a  given  tolerance  above  the  target  apogee.  Second,  the 
predicted  final  energy  indicated  that  the  vehicle  would  be  on  an  elliptical,  not  hyperbolic, 
trajectory.  So  as  to  incorporate  generic  features,  modeling  of  a  separate  lift  down  phase  was 
replaced  by  using  a  lift-down  condition  to  determine  phase  change  and  biasing  to  the  same  lift- 
down  condition  during  targeting.  Also,  use  of  a  heuristic  sensitivity  to  calculate  the  first 
corrected  bank  angle  was  replaced  by  a  simple  ‘smart  guessing’  algorithm.  Finally,  heuristic 
lateral  corridor  boundaries  were  replaced  by  boundaries  based  on  percentage  of  forward  velocity. 

Analysis  of  the  resulting  entry  corridor  revealed  that  the  enhanced  algorithm  generated  a 
mean  improvement  of  261%  over  the  PredGuid  corridor.  Most  of  the  gain  occurred  for  steeper 
flight  path  angles.  Results  also  indicated  that  the  enhanced  algorithm  yielded  lower  maximum 
accelerations  and  comparable  heating  rates,  heating  loads,  and  AV  to  raise  perigee.  These  results 
are  intended  to  provide  a  starting  point  for  further  enhancement  and  applicability  to  interplanetary 
travel. 
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Chapter  1 
Introduction 

1.1  Aerocapture 

Aerocapture  is  an  expeditious,  fuel-efficient  method  of  orbit  insertion.  Payload  and 
vehicle  mass  are  often  limited  by  the  amount  of  chemical  propellant  required  for 
propulsive  orbit  insertions.  As  a  result,  methods  reducing  the  required  propellant  are 
highly  desirable.  The  most  readily  available,  non-propulsive,  means  of  slowing  a  vehicle 
from  hyperbolic  speeds  to  low  orbital  speeds  is  atmospheric  drag.  Variants  of  this 
aeroassist  approach  have  surfaced  over  the  last  several  decades.  One  popular  variant, 
aerobraking,  uses  a  propulsive  maneuver  to  insert  the  vehicle  into  a  high  elliptical  orbit. 
Multiple  passes  through  the  atmosphere  slowly  remove  velocity  in  small  increments  until 
the  orbit  is  circular.  While  aerobraking  can  considerably  reduce  the  amount  of  propellant 
required,  a  significant  amount  is  still  necessary  for  the  initial  orbit  insertion. 

Aerocapture  completely  eliminates  the  need  for  propulsive  insertion  by  use  of  a  single 
pass  through  the  atmosphere  to  slow  a  vehicle  from  the  high-energy  approach  trajectory 
to  a  low-energy  orbit.  Drag  forces  remove  energy  until  the  vehicle  is  slowed  below 
escape  velocity.  Unlike  aerobraking,  there  is  no  intermediary  elliptical  orbit.  Thruster 
firings  are  used  in  the  process,  but  they  are  limited  to  controlling  orientation  and  provide 
no  reduction  in  velocity.  Figure  1.1  highlights  differences  between  the  two  techniques. 


Aerobraking  Aerocapture 


Figure  1.1:  Aerocapture  vs.  Aerobraking 
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Numerous  studies  have  been  conducted  on  the  benefits  of  aerocapture.  Reference  [4] 
compared  aerobraking  and  aerocapture  for  the  Mars  Orbiter.  These  studies  indicated  a 
reduction  in  total  vehicle  mass  from  1600  kg  to  1 100  kg  and  a  reduction  in  time  required 
for  orbit  insertion  from  3  months  to  3  hours  through  use  of  aerocapture. 

Even  with  these  promising  figures  and  availability  of  the  technology  since  the  Apollo 
missions,  aerocapture  has  never  been  used.  The  Apollo  Program  had  a  technique  on 
paper  similar  to  aerocapture,  allowing  entry  capsules  to  enter  the  atmosphere  and  exit  at 
near-orbital  velocity  as  a  method  of  landing  area  weather  avoidance.  However,  weather 
avoidance  was  never  required  during  the  program,  so  the  technology  was  developed  and 
man-rated  but  not  demonstrated  [3]. 

1.2  Previous  Research 

Several  proposals  for  an  aerocapture  flight  demonstration  have  emerged  in  the  last 
decade,  including  the  Aeroassist  Flight  Experiment  in  the  early  1990’s  and  more  recently, 
the  joint  NASA-CNES  Mars  Sample  Return  Mission  [3].  These  proposals  led  to  the 
development  of  several  important  algorithms  used  in  simulations  to  verify  aerocapture. 
They  are  an  analytic  predictor-corrector,  a  numeric  predictor-corrector,  a  terminal  point 
controller,  and  an  energy  controller. 

1.2.1  Analytic  Predictor-Corrector 

The  analytic  predictor-corrector  is  the  second  phase  of  a  two-phase  aerocapture 
maneuver.  Dividing  the  aerocapture  maneuver  into  two  phase  allows  separate  control  of 
trajectory  loads  and  apogee  targeting.  An  equilibrium  glide  phase  comprises  the  first  part 
of  the  trajectory  and  the  analytic  predictor-corrector  is  used  for  the  exit  phase.  In  the 
equilibrium  glide  phase,  the  commanded  bank  angle  is  modeled  with  a  linear  second- 
order  differential  equation  for  altitude.  After  the  vehicle  has  been  slowed  to  a  particular 
velocity,  it  transitions  to  the  exit  phase. 

The  analytic  integration  used  by  the  predictor-corrector  reduces  the  onboard  computer 
requirements  as  compared  to  a  numerical  method.  Analytic  integration  is  possible  by 
using  altitude  rate,  h,  as  the  only  control  variable.  Reference  [1]  contains  a  detailed 
derivation  of  the  analytic  equations  of  motion.  These  analytic  equations  predict  the 
vehicle  relative  velocity  at  exit,  assuming  a  constant  altitude  rate.  Relative  velocity  is 
converted  to  inertial  velocity  and  the  predicted  altitude  rate  and  inertial  velocity  at  exit 
are  used  to  predict  the  resulting  apogee.  The  algorithm  iterates  on  the  initial  altitude  rate 


18 


until  the  desired  apogee  is  achieved.  This  altitude  rate  is  used  to  calculate  commanded 
bank  angle  with  an  equation  similar  to  the  control  equation  from  the  equilibrium  glide 
phase. 

Lateral  control  determines  the  sign  of  the  bank  angle  in  order  to  keep  the  orbital  plane  as 
close  as  possible  to  the  target  orbital  plane.  Since  there  is  only  one  control  variable,  bank 
angle,  it  is  impossible  to  simultaneously  null  both  position  and  velocity  errors.  However, 
it  is  possible  to  minimize  both  errors  by  nulling  the  wedge  angle,  5,  which  is  the  angle 
between  the  actual  and  desired  orbital  planes. 

References  [3]  and  [2]  provide  a  more  detailed  description  of  this  analytic  predictor- 
corrector  method. 

1.2.2  Numeric  Predictor-Corrector 

The  numeric  predictor-corrector  algorithm  controls  the  orientation  of  the  lift  vector  about 
the  relative  velocity  vector  by  altering  the  bank  angle.  The  more  deceleration  required, 
the  deeper  the  vehicle  will  penetrate  into  the  atmosphere.  The  algorithm  numerically 
integrates  the  current  position  and  velocity  vectors  forward  to  atmospheric  exit  using  a 
constant  bank  angle  assumption.  At  each  integration  step,  gravitational  and  aerodynamic 
accelerations  are  calculated  using  simplified  models.  The  predicted  states  are  used  to 
calculate  the  resulting  apogee.  Bank  angle  is  then  adjusted  to  null  the  target  apogee  miss 
and  this  cycle  is  continued  until  the  final  apogee  is  within  acceptable  tolerances. 

Lateral  control,  as  with  the  analytic  predictor-corrector,  determines  the  sign  of  the  bank 
angle  in  order  to  keep  the  orbital  plane  as  close  as  possible  to  the  target  orbital  plane.  As 
discussed  in  the  last  section,  it  is  impossible  to  simultaneously  null  both  position  and 
velocity  errors.  This  algorithm  controls  the  velocity  error  due  to  the  short  duration  of  the 
aerocapture  trajectory. 

The  numeric  predictor-corrector  is  described  in  reference  [1]  and  is  the  basis  for  the 
aerocapture  algorithm  used  in  this  thesis. 

1.2.3  Terminal  Point  Controller 

The  terminal  point  controller  algorithm  attempts  to  drive  the  vehicle  to  follow  a 
predetermined  reference  trajectory  to  a  fixed  terminal  point  or  set  of  terminal  conditions. 
The  most  crucial  part  of  using  a  terminal  point  controller  is  generating  an  optimal 
reference  trajectory.  Simple  feedback  guidance  schemes  correct  for  dispersions  and  other 
errors  along  the  trajectory. 


19 


As  with  the  numeric  and  analytic  predictor-corrector  methods,  the  terminal  point 
controller  maintains  lateral  control  through  the  use  of  predetermined  limits  on  corridor 
error. 

The  terminal  point  controller  is  described  in  references  [7]  and  [8]. 

1.2.4  Energy  Controller 

In  the  energy  controller  algorithm,  the  vehicle  energy  is  controlled  to  a  targeted  energy 
state  (determined  by  the  target  apogee)  by  altering  the  energy  gain.  The  energy  gain  is 
calculated  by  taking  the  ratio  of  energy  rate  (a  function  of  drag)  to  energy  error.  The  gain 
is  controlled  so  that  the  error  and  the  energy  rate  approach  zero.  Energy  gain  is  translated 
into  altitude  rate  and  the  bank  angle  is  determined  using  an  analytical  equation  for  vehicle 
vertical  acceleration. 

Lateral  control  is  similar  to  that  used  in  the  other  algorithms.  If  the  out  of  plane  velocity 
exceeds  a  given  deadband,  the  vehicle  executes  a  roll  reversal  to  reduce  the  out  of  plane 
velocity  error. 

This  algorithm  is  described  in  detail,  including  all  relevant  derivations,  in  reference  [1]. 


1.3  Problem  Definition  and  Thesis  Objective 

Over  the  last  two  decades,  The  Charles  Stark  Draper  Laboratory  (CSDL)  has  investigated 
the  aerocapture  concept.  In  1984,  John  Higgins  developed  the  original  numeric 
predictor-corrector  targeting  guidance  algorithm  for  Earth  orbit  transfer  applications  [1]. 
This  algorithm  has  been  periodically  updated  for  use  in  various  aerocapture  guidance 
proposals  but  has  its  limitations.  An  aerocapture  vehicle  is  capable  of  capturing  within  a 
certain  range  of  flight  path  angles.  The  current  application  of  Higgins’  algorithm, 
PredGuid,  is  very  useful  for  low-energy,  shallow  trajectories.  As  the  energy  of  the  orbits 
increases  and  the  entry  flight  path  angles  become  steeper,  the  portion  of  the  theoretical 
corridor  successfully  captured  by  the  algorithm  decreases.  The  theoretical  corridor  is 
defined  by  the  full  lift  up  and  full  lift  down  flight  path  angle  boundaries  as  well  as 
structural  and  thermal  constraints.  The  capturable  corridor  is  the  percentage  of  the 
theoretical  corridor  for  which  the  algorithm  allows  the  vehicle  to  aerocapture  successfully 
within  target,  structural,  and  thermal  constraints.  These  corridors  will  be  further 
described  in  Chapter  4. 
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One  major  reason  for  the  breakdown  of  Higgins’  algorithm  at  high-energy  cases  is  the 
inability  of  the  algorithm  logic  to  handle  hyperbolic  trajectories.  This  algorithm 
calculates  updated  bank  angle  guesses  based  on  a  slope  calculated  from  previous  guesses. 
Hyperbolic  trajectories  have  negative  values  for  apogee  and  semi-major  axis  which  are 
incorrectly  interpreted  by  the  algorithm  as  Tow’  misses.  The  second  reason  for  high- 
energy  case  failure  is  the  constant  bank  angle  assumption.  The  algorithm’s  only  recourse 
for  reducing  energy  is  digging  deeper  into  the  atmosphere  which  in  turn  raises 
acceleration  loads  to  unacceptable  levels. 

In  his  1988  Master  of  Science  thesis,  Doug  Fuhry  added  an  energy  management 
capability  to  Higgins’  algorithm  with  a  numerical  predictor-corrector  entry  phase 
targeting  into  a  constant  altitude  cruise  phase  [6].  At  a  specified  velocity  the  vehicle 
would  then  transition  to  Higgins’  targeting  algorithm.  These  improvements  allowed 
enhanced  aerocapture  coverage  for  applicability  at  Mars.  However,  the  algorithm  was 
fairly  complex,  containing  three  phases  for  the  aerocapture  maneuver. 

Although  Higgins’  numeric  predictor-corrector  method  is  generic  in  nature  and  can  be 
used  for  any  planet  or  moon  with  an  atmosphere,  the  algorithm  was  designed  for  use  at 
Earth.  Similarly,  Fuhry’ s  method,  like  those  discussed  in  the  previous  section,  while 
generic  in  concept,  was  specialized  for  aerocapture  at  Mars.  Through  empirical 
observations  from  thousands  of  test  cases,  these  algorithms  incorporated  new  heuristic 
approximations  and  features  for  optimization.  Use  of  the  existing  algorithms  for  a 
different  case  would  involve  empirically  determining  and  I-loading  a  large  number  of 
factors.  However,  a  more  generic  algorithm,  capable  of  expanding  the  capturable 
corridor  and  easily  applicable  to  different  planets  and  moons  with  atmospheres,  is 
desirable. 

This  thesis  seeks  to  enhance  the  numerical  predictor-corrector  aerocapture  guidance 
algorithm  (PredGuid)  by  implementing  a  single  energy  management  phase  prior  to 
targeting,  developing  a  generic  method  of  transitioning  between  the  energy  management 
and  targeting  phases,  and  replacing  other  heuristic  features  with  more  generic  features. 
The  resulting  flight  path  angle  entry  corridor  will  be  compared  to  the  flight  path  angle 
entry  corridor  of  the  original  algorithm.  Additionally,  comparisons  will  be  made  for 
various  characteristics  of  the  trajectories  including  maximum  g  load,  heating  rate,  heating 
load,  and  AV  to  raise  perigee.  The  resulting  algorithm  is  intended  to  provide  a  starting 
point  for  further  enhancement  for  applicability  to  interplanetary  travel. 
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1.4  Thesis  Overview 


This  thesis  contains  details  and  results  of  the  design  and  testing  of  the  enhanced 
aerocapture  guidance  algorithm  discussed  above. 


Chapter  2,  Equations  of  Motion,  describes  the  equations  used  in  the  simulation  dynamics. 
Reference  frames,  coordinate  systems,  assumptions,  environmental  models,  and  vehicle 
models  used  in  the  dynamics  are  described. 


Chapter  3,  Guidance  Design,  describes  the  algorithm  used  in  the  aerocapture  guidance. 
Design  factors,  guidance  phases,  phase  change  criteria,  and  important  alterations  to  the 
original  algorithm  are  described. 


Chapter  4,  Capture  Envelope  and  Test  Case  Determination,  describes  the  theoretical  and 
capturable  corridors.  Corridor  approximations,  vehicle  constraint  concerns,  and 
determination  of  initial  conditions  are  described. 


Chapter  5,  Results,  discusses  the  results  of  the  various  cases  tested.  Comparisons  to  the 
original  algorithm  as  well  as  evaluation  of  the  enhanced  algorithm  are  described. 

Chapter  6,  Conclusions,  discusses  conclusions  of  this  study  as  well  as  suggestions  for 
future  study. 
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Chapter  2 

Equations  of  Motion 

2.1  Overview 

Equations  of  motion  model  the  environment  to  provide  an  accurate  prediction  of  the 
vehicle  trajectory  within  the  simulation.  Components  of  the  equations  of  motion  are 
computed  in  different  reference  coordinate  frames,  making  it  necessary  to  define  both  the 
frames  and  relevant  rotation  and  transformation  matrices  between  them.  This  chapter 
provides  a  brief  guide  to  the  reference  frames  and  rotations  used  in  this  thesis  as  well  as 
the  basic  equations  of  motion  used  in  the  environment  models. 

2.2  Reference  Coordinate  Frames 

Inertial  Reference  Frame  ( ii,ji,ki  ):  a  non-rotating  Earth-centered  coordinate  system 
in  which  the  origin  lies  at  the  center  of  the  Earth.  The  it  axis  points  through  zero 
longitude  at  time  zero,  the  k,  axis  points  through  the  North  Pole,  and  the  ji  axis 
completes  the  right-handed  coordinate  system. 


A  A  A  %  §  # 

Local  Horizontal  Reference  Frame  ( ih  ,jh,kh ):  a  coordinate  system  with  the  origin  at 
the  vehicle  center  of  gravity  (CG).  The  jh  axis  is  in  the  direction  of  the  angular 
momentum  vector  calculated  from  the  relative  velocity  (Vrel  xR).  The  kh  axis  points 
towards  the  center  of  the  Earth  along  the  vehicle  inertial  position  vector,  and  ih 
completes  the  right-handed  coordinate  system. 


Velocity  Reference  Frame  ( iv,jv,kv  ):  a  coordinate  system  with  the  origin  at  the  vehicle 
CG  and  the  iv  axis  pointing  along  the  vehicle  relative  velocity  vector.  The  jv  axis 
remains  in  the  local  horizontal  and  the  kv  axis  completes  the  right-handed  coordinate 
system. 
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Body  Reference  Frame  (ib,jb,kb  ):  a  coordinate  system  with  the  origin  at  the  vehicle 
CG.  The  ib  axis  points  along  the  longitudinal  axis  of  the  vehicle.  The  jb  axis  is  positive 
out  the  ‘right  wing’.  The  kb  axis  completes  the  right-handed  coordinate  system  and 
points  positive  downward  (towards  Earth). 


Stability  Reference  Frame  {is,js,ks  ):  a  coordinate  system  with  the  origin  at  the  vehicle 
CG  with  the  is  axis  along  the  projection  of  the  velocity  vector  (V)  onto  the  body  ib-kb 
plane.  Assuming  zero  sideslip,  the  is  and  iv  axes  are  coincident.  The  js  axis  is 
coincident  with  the  jb  axis  and  the  ks  axis  completes  the  right-handed  coordinate 
system. 

2.3  Coordinate  Transformations 


Vector  transformations  are  accomplished  by  taking  the  dot  product  of  a  vector  in  a  given 
frame  with  the  basis  vectors  of  a  desired  reference  frame.  This  operation  can  be  written 
in  matrix  form  as  a  multiplication  of  a  transformation  matrix  (denoted  Ta2b)  with  a  vector. 
The  transformation  matrix  contains  the  unit  vectors  of  the  given  frame  expressed  in  the 
desired  reference  frame.  For  example,  a  transformation  matrix  from  reference  frame  A  to 
reference  frame  B  would  contain  the  basis  vectors  of  A  written  in  the  B  frame  as  shown 


in  Eq.  2.1. 
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Reference  frames  are  generally  related  by  an  easily  identifiable  angle  rotated  about  one 
axis.  As  a  result,  transformation  matrices  are  most  frequently  seen  as  combinations  of 
sines  and  cosines  on  two  axes  while  the  third  axis  remains  aligned  with  the  given  frame. 
Such  transformation  matrices  are  referred  to  as  Euler  rotations  and  are  extremely  useful 
in  aerospace  applications.  Roll,  pitch,  and  yaw  Euler  rotations  are  the  most  common. 
Rotation  or  transformation  in  the  opposite  direction  is  accomplished  by  simply  using  the 
transpose  of  the  relevant  matrix.  Eq.  2.2  and  2.3  demonstrate  this  principle. 
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(2.3) 


Velocity  to  Inertial  Transformation:  The  transformation  matrix  between  the  velocity 
and  inertial  reference  frame  contains  the  unit  vectors  of  the  velocity  reference  frame 
written  in  the  inertial  frame. 
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where: 
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Bank  Angle  Rotation:  The  stability  reference  frame  is  the  velocity  reference  frame 
rotated  by  only  the  bank  angle  <()  about  their  common  i  axes.  This  relationship  is  shown 
in  Figure  2.1  below  and  the  transformation  matrix  between  the  frames  is  given  in  Eq.  2.8. 
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Figure  2.1:  Bank  Angle  Rotation 
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Angle  of  Attack  Rotation:  The  body  reference  frame  is  the  stability  reference  frame 
rotated  by  the  angle  of  attack  a  about  their  common  j  axes.  This  relationship  is  shown 
in  Figure  2.2  below  and  the  transformation  matrix  between  the  frames  is  given  in  Eq.  2.9. 


Figure  2.2:  Angle  of  Attack  Rotation 
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Aerodynamic  accelerations  act  on  the  body  in  the  stability  frame.  However,  propagation 
of  the  states  is  accomplished  in  the  inertial  frame.  Therefore,  a  single  transformation 
matrix  to  rotate  between  the  stability  and  inertial  reference  frames  is  desirable.  This 
transformation  matrix  can  be  expressed  as  a  product  of  the  above  matrices  and  is  given  in 
Eq.  2.10. 

T  =T  T  (2.10) 

*  s2i  1  vli1  slv  ^  ^ 


Similarly,  the  body  to  inertial  transformation  can  be  expressed  as: 


T  _  T  T  T 
1b2i  1v2i1s2v1b2s 


(2.11) 


2.4  Environment  Models 


Environment  models  describe  the  surroundings  to  which  the  vehicle  is  exposed.  In  this 
thesis,  there  are  two  basic  models  that  must  be  considered:  atmosphere  and  gravity. 
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2.4.1  Earth  Atmosphere  Models 


Atmosphere  models  are  used  to  calculate  the  density  of  the  atmosphere  at  a  given  height 
above  the  surface.  This  density  is  important  for  lift  and  drag  calculations.  This  thesis 
uses  the  US  Standard  Atmosphere,  1976  for  the  dynamic  simulation  [10]. 


2.4.2  Earth  Gravity  Model 


Gravity  models  are  used  to  calculate  the  acceleration  due  to  gravity  of  the  vehicle  at  a 
given  position  with  respect  to  the  planet.  Earth-based  gravity  models  can  be  simple  in 
nature  and  increase  in  complexity  depending  on  the  desired  accuracy.  For  this  thesis, 
normal  conical  acceleration  and  J2  effects  were  considered: 
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Eq.  2.12  is  derived  in  reference  [9]. 


2.5  Equations  of  Motion 


The  full  equations  of  motion  for  an  aerocapture  vehicle  are  derived  using  Newton’s 
Laws.  Several  assumptions  are  used  to  simplify  the  derivation.  These  assumptions  are  as 
follows: 

■  The  vehicle  is  symmetrical. 

■  The  vehicle  has  constant  mass  (i.e.  fuel  bum  is  not  accounted  for). 

■  The  vehicle  produces  no  thrust. 

■  There  are  no  aerodynamic  moments  (vehicle  is  statically  trimmed). 

■  There  are  no  side  forces  (no  sideslip,  (3). 

With  these  assumptions,  acceleration  acting  on  the  vehicle  is  comprised  of  two  parts: 
acceleration  due  to  gravity  and  acceleration  due  to  aerodynamic  forces: 

a  =  ag+aa  (2.13) 

Aerodynamic  acceleration  is  further  broken  down  into  two  components:  acceleration  due 
to  drag  and  acceleration  due  to  lift. 

aB=ad+a,  (2.14) 
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Since  this  simulation  was  developed  without  a  solid  vehicle  description,  general 
parameters  such  as  ballistic  coefficient  (CB)  and  lift-to-drag  ratio  (L/D) were  sufficient 
to  calculate  the  aerodynamic  accelerations.  The  force  due  to  drag  ( FD )  acting  on  a 
vehicle  is  defined  using  the  vehicle  coefficient  of  drag  ( Cd ),  dynamic  pressure  {(/ ) ,  and 
planform  (or  reference)  area  (S ) . 

FD=CDqS  (2.15) 


(2.16) 


Eq.  2.17  defines  ballistic  coefficient: 


CB  = 


m 

~C^S 


(2.17) 


Using  these  definitions  and  Newton’s  Law  of  gravitation  with  constant  mass. 


Fd  =  maD 


acceleration  due  to  drag  ( aD )  can  be  derived  as: 
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Then  acceleration  due  to  lift  (a,  )  is  found  simply: 
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ai  =ClD - 

1  D 


(2.20) 


Drag  acceleration  acts  in  the  direction  opposite  the  relative  velocity  while  lift  acceleration 
acts  perpendicular  to  relative  velocity.  Written  in  the  stability  frame,  the  resulting 
aerodynamic  acceleration  is: 
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For  integration,  acceleration  must  be  rotated  into  the  inertial  frame. 
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Yielding: 

a0i  =  ~aDiVi  +aL  sin(^>V2  -aL  cos (<j>\ 

aa.  =  -aDjv>  +  aL  sin -  a,  cos {<j))j^  (2.23) 
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Acceleration  due  to  gravity  can  also  be  expressed  as  a  sum  of  components:  conical 
acceleration  and  acceleration  due  to  a  nonspherical  body. 

ag=ac+ans  (2.24) 


Conical  acceleration  is  expressed  as: 
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(2.25) 


The  nonspherical  terms  can  be  added  to  the  conical  equation  to  account  for  the  Earth’s 
oblateness,  non-uniform  density,  etc.  For  the  purposes  of  this  thesis  the  J2  terms  are 
sufficient  and  the  resulting  total  gravitational  acceleration  is  the  expression  from  Eq.  2.12 
above.  Note  that  this  acceleration  is  already  written  in  the  inertial  frame  and  no  further 
rotation  is  required. 


Finally,  the  total  acceleration  in  the  inertial  frame: 
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Integrating  acceleration  once  yields  the  change  in  velocity  in  the  inertial  frame. 
Integrating  twice  yields  the  change  in  position  in  the  inertial  frame.  Appendix  B  contains 
a  detailed  verification  that  these  equations  of  motion  were  properly  implemented  into  the 
system  dynamics. 


2.6  Vehicle  Models 


The  vehicle  used  in  this  study  was  described  solely  by  mass,  planform  area,  lift-to-drag 
ratio,  and  coefficient  of  drag.  The  values  used  in  this  thesis  are  representative  of  a 
conceivable  aerocapture  vehicle. 

2.7  Vehicle  Properties 

The  following  values  were  used  for  mass,  planform  area,  and  drag  coefficient: 
m  =  15.4783  slugs  (225.89  kg) 

S  =  12.163  ft2  (l,13  m2) 
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Cd  =1.4286 

These  values  yield  a  ballistic  coefficient,  CB  =  0.891157  slugs  /ft2  (l39.93  kg/m  ). 

The  algorithm  was  highly  sensitive  to  the  lift-to-drag  ratio.  As  a  result,  for  testing,  lift-to- 
drag  ratio  was  varied  between  0.25  and  1.5  in  order  to  gain  a  better  understanding  of  the 
range  of  possible  trajectories. 

2.7.1  Aerodynamic  Heating 


Heating  rate  and  heating  load  are  important  factors  in  aerocapture  vehicle  performance. 
The  convective  heating  rate  is  calculated  empirically  from  an  equation  in  reference  [5]. 
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where  Rn  is  the  vehicle  nose  radius,  pSL  is  the  density  at  the  surface  of  the  earth,  and  ve 
is  the  reference  spherical  velocity  given  by: 

(2.28) 


The  following  values  were  used: 

Rn  =  1.96764  ft  (0.60  m) 

ve  =25936.241/// 5  (7.91  km/s) 

pSL  =  0.00237688  slugs/  ft 3  (l.22  kg/m 3) 

Nose  radius  was  calculated  by  assuming  that  the  planform  area  of  the  vehicle  is  a  circle 
and  the  radius  of  that  circle  is  the  nose  radius. 

The  heat  load  due  to  aerodynamic  heating  is  found  by  integrating  the  heating  rate  over 
time: 

Q(t)=  [Q(r}ir  (2.29) 
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Chapter  3 
Guidance  Design 

3.1  Overview 

The  starting  point  for  this  thesis  was  a  targeting  algorithm,  known  as  PredGuid,  based  on 
Higgins’  numeric  predictor-corrector  [1],  The  stated  objective  of  this  thesis  was  to 
enhance  this  existing  targeting  algorithm  by  implementing  an  energy  management  phase 
prior  to  targeting,  developing  a  generic  method  of  transitioning  between  the  energy 
management  and  targeting  phases,  and  replacing  other  heuristic  features  with  more 
generic  features.  A  detailed  description  of  PredGuid,  with  all  the  mission-specific 
heuristics,  is  contained  in  Appendix  A.  This  chapter  describes  the  enhanced  guidance 
algorithm  and  phase  change  logic  in  entirety  and  specifically  details  the  enhancements  to 
Higgins’  targeting  algorithm. 

The  enhanced  algorithm  flow  (depicted  in  Figure  3.1)  is  identical  to  PredGuid.  However, 
in  the  enhanced  algorithm,  Bank  Angle  Determination,  which  was  originally  only  a 
targeting  phase,  is  now  further  broken  down  into  two  phases,  Energy  Management  and 
Targeting,  discussed  in  sections  3.2  and  3.3,  respectively.  Lateral  Control,  discussed  in 
section  3.4,  and  Command  Incorporation,  discussed  in  section  3.5  retain  their  PredGuid 
objectives,  but  have  been  modified  appropriately  to  remove  heuristics  and  support  the 
addition  of  the  Energy  Management  Phase.  Initialization  and  Aerodynamic  Properties 
were  not  substantially  modified  from  PredGuid.  Changes  were  limited  to  cosmetics  and 
substitution  of  the  US  Standard  Atmosphere,  1976  in  place  of  the  original  exponential 
atmospheric  model.  Initialization  is  executed  once  on  the  first  guidance  cycle. 
Aerodynamic  Properties  calculates  estimated  density,  lift-to-drag  ratio,  and  drag 
coefficient.  These  values  are  updated  throughout  the  sensible  atmosphere,  and  for  a  short 
time  before  guidance  begins  cycling,  to  ensure  that  guidance  is  receiving  good  data. 
Appendix  A  contains  a  detailed  description  of  these  procedures. 
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Figure  3.1:  Guidance  Algorithm  Flow 

During  aerocapture,  the  vehicle  steers  solely  by  rotating  the  lift  vector  about  the  relative 
velocity  vector.  Guidance  assumes  the  vehicle  makes  coordinated  turns  and  trims  to  a 
constant  angle  of  attack.  Therefore,  bank  angle  (<f)  is  the  only  control  factor.  Positive 
bank  is  defined  as  bank  to  the  right  according  to  the  right  hand  rule.  Bank  angle  can  vary 
between  ±180°  where  0°is  full  lift  up  (away  from  Earth)  and  180°  is  full  lift  down 
(towards  Earth).  However,  bank  control  is  typically  limited  between  15°  and  165°  (in 
either  direction)  by  the  guidance  algorithm  in  order  to  constantly  maintain  some  control 
authority  for  out-of-plane  corrections.  The  algorithm  assumes  the  in-plane  and  out-of¬ 
plane  velocity  controls  are  decoupled.  Therefore,  the  bank  angle  required  to  reach  the 
target  apogee  is  computed  independently  of  the  bank  angle  direction  to  remain  in  the 
desired  orbital  plane. 

Guidance  is  initialized  at  entry  interface  (El)  by  running  the  energy  management 
guidance  cycle  once.  Guidance  does  not  cycle  again  until  the  vehicle  is  in  the  sensible 
atmosphere  in  order  to  prevent  control  corrections  when  there  is  insufficient  atmospheric 
density  to  have  reasonable  command  authority.  Entry  into  the  sensible  atmosphere  is 
determined  by  aerodynamic  acceleration,  calculated  in  Eq.  3.1: 

K<™|  (3.1) 

u aero  v  7 

s © 

While  the  vehicle  is  experiencing  loads  above  the  acceleration  minimum  (0.075  g’s), 
guidance  cycles  at  a  frequency  of  0.5  Hz. 

3.2  Energy  Management  Phase 

The  purpose  of  the  energy  management  phase  is  to  deplete  sufficient  energy  to  allow 
targeting  of  the  exit  conditions  in  the  next  phase.  Several  approaches  have  been  explored 
including  maintaining  constant  altitude,  specified  dynamic  pressure,  or  reference  drag. 
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All  of  these  approaches  adhered  to  a  similar  format,  containing  an  altitude  rate  damper 
and  a  reference  following  term.  The  algorithm  used  in  this  thesis  attempts  to  null  out  the 
altitude  rate  while  maintaining  a  reference  drag.  This  method  was  chosen  for  ease  of 
implementation.  The  reference  drag  term  can  easily  be  replaced  with  a  reference  altitude 
term  or  dynamic  pressure  term  and  achieve  similar  results.  While  maintaining  the 
constant  cruise  condition,  the  guidance  runs  periodic  checks  to  determine  when  to  initiate 
the  targeting  phase.  The  two  major  components  of  the  Energy  Management  Phase  are 
discussed  in  the  following  sections.  The  Energy  Management  Phase  is  depicted  in  Figure 
3.2. 


Figure  3.2:  Energy  Management  Phase 


3.2.1  Constant  Condition  Cruise 


The  constant  condition  cruise  phase  modeled  in  this  thesis  draws  upon  a  conglomeration 
of  the  algorithms  presented  in  references  [2],  [1]  and  [6],  The  vehicle  response  can  be 
modeled  as  a  second-order  spring/mass/damper  system.  Beginning  with  Newton’s 
Second  Law,  the  acceleration  of  the  vehicle  in  the  inertial  reference  frame  can  be  written: 


F  _  d2r 
m  dt 1 


(3.2) 


Aerodynamic  forces  are  most  often  expressed  in  local  coordinate  systems.  Because  the 
local  reference  frame  is  rotating,  we  must  include  the  coriolis  terms  in  the  derivatives. 
The  first  derivative  of  position  in  the  inertial  frame  is  then: 

dr  _  dr 
dt  j  dt 


+'cohy.r  (3.3) 

h 


where  h  denotes  the  local  horizontal  reference  frame  from  Chapter  2.  Taking  the  second 
derivative, 
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+'3h  x(,a>h  xr) 


(3.4) 
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Position  and  relative  velocity  in  the  inertial  frame  are  defined  in  the  local  horizontal 
frame  in  Eq.  3.5  And  3.6,  respectively: 


r1  =  -rhk. 


v',  -  vrcl  cos {y\  -  vrd  sin {y)kh 


(3.5) 

(3.6) 


where  y  is  the  flight  path  angle.  Angular  velocity  of  the  local  horizontal  frame  with 
respect  to  the  inertial  frame  is  defined  in  Eq.  3.7. 
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For  small  y, 


(3.8) 

(3.9) 


For  the  energy  management  phase,  we  are  only  concerned  with  forces  and  acceleration  in 
the  kh  direction.  Substituting  the  simplifications  above  into  Eq.  3.2  and  Eq.  3.4,  we  find 
that,  in  the  kh  direction, 


F_ 

m 


(3.10) 


Lift  and  gravity  are  the  only  forces  in  the  kh  direction,  therefore, 

F  =  -Lcos W  +  g  =  _..  +  ^,  (3.11) 

mm  r 


Using  the  definitions  of  drag  and  ballistic  coefficient, 

D  =  CDqS 


(3.12) 


C„  = 


m 

C^S 


acceleration  in  the  kh  is  presented  in  Eq.  3.14. 
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Setting  this  acceleration  equal  to  zero,  the  cosine  of  bank  required  to  achieve  zero  altitude 
rate  is: 


cos 


<f>\ 


1 

r=0  q 

yDy 

(3.15) 


The  commanded  bank  angle  from  guidance  is  calculated  by  taking  the  value  of  bank 
required  for  zero  altitude  rate  and  adding  an  altitude  rate  damper  and  reference  drag 
following  terms: 


^  A  tr  r  ,  ^  l D~Dref ) 

cos  =  cos  <f>\  -Kr—  +  KD - = - 

r~u  q  q 


(3.16) 


Substituting  Eq.  3.14  into  Eq.  3.16  for  the  cosine  term,  the  second-order  equation 
becomes: 

r  +  ——Kfr  +——KD(D-Dref)=  0  (3.17) 

CBD  r  CbD  dV  refJ 


The  following  equation  for  reference  drag  (with  a  gain  value  of  K  =  1.5)  was  taken  from 
reference  [2], 


Dref  =  K 


V 

— g 


Cr 


(3.18) 


The  general  equation  of  a  second  order  spring/mass/damper  system  is: 

x  +  2  o)n  ^pc  +  (o2n  —  0 


(3.19) 


Applying  this  general  equation  to  Eq.  3.17,  natural  frequency  and  damping  ratio  are: 

co2=— — Kd  (3.20) 

"  CbD  d 

2co£  =  ——K,  (3.21) 

"  CB  D  r 

The  following  values  were  extracted  from  references  [1]  and  [6]  and  worked  well  for  this 
application: 

con  =  0.06  rad  /  s 
£  =  1.5 

Note  that  since  the  lift-to-drag  ratios  were  varied  in  this  study,  the  resulting  gain  values 
are  functions  of  L/D. 
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Kd  =  f(L/D) 

K,=f{L/D) 

3.2.2  Phase  Change  Conditions 

The  vehicle  remains  in  the  energy  management  phase  until  sufficient  energy  has  been 
depleted  to  allow  targeting.  This  condition  is  quantified  by  determining  if,  at  the  given 
point  in  the  trajectory,  the  vehicle  could  fly  a  constant  bank  angle  of  110°  for  the 
remainder  of  the  trajectory  and  have  the  resulting  apogee  be  below  or  within  a  given 
tolerance  above  the  target  apogee.  The  bank  angle  of  110°  was  not  chosen  arbitrarily. 
Extensive  research  completed  in  past  studies  indicates  that  flying  a  slightly  lift  down 
bank  angle  when  exiting  the  atmosphere  raises  the  resulting  periapse  altitude  [1]. 
Reducing  the  amount  of  AV  required  to  raise  the  periapse  altitude  after  aerocapture 
greatly  reduces  the  fuel  budget. 

One  more  condition  is  required  in  order  to  ensure  guidance  does  not  change  phases 
prematurely.  Hyperbolic  orbits  have  negative  values  of  semi -major  axis.  If  the  predicted 
final  position  and  velocity  are  still  on  a  hyperbolic  orbit,  the  radius  of  apogee  will  be 
calculated  to  be  a  negative  number.  Obviously,  a  hyperbolic  trajectory  is  not  a  low 
trajectory;  however,  a  negative  apogee  will  be  interpreted  by  targeting  as  low.  If 
guidance  exits  energy  management  too  early  based  on  the  Tow’  predicted  apogee,  the 
orbit  will  remain  hyperbolic  and  the  targeting  algorithm  will  be  unable  to  resolve  this 
condition  to  reach  the  target  apogee.  To  avoid  premature  exit,  guidance  checks  the 
energy  of  the  predicted  exit  conditions.  Energy  will  be  positive  on  a  hyperbolic  orbit, 
zero  on  a  parabolic  orbit,  and  negative  on  an  elliptical  orbit.  The  algorithm  does  not 
change  phases  until  the  energy  of  the  predicted  exit  conditions  is  negative.  This  ensures 
that  the  initial  phase  check  does  not  give  a  false  indication  of  requiring  a  phase  change. 

s  =  -v2  -tL  =  -IL  (3.22) 

2  r  2a 

The  phase  change  check  executes  at  a  slower  frequency  than  the  main  guidance,  running 
at  a  frequency  of  0.1667  Hz.  Figure  3.3  depicts  the  functional  flow  of  the  phase  change 
logic. 
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Figure  3.3:  Phase  Change  Functional  Flow 


3.3  Targeting  Phase 

The  targeting  portion  of  the  guidance  algorithm  relies  on  a  numeric  predictor-corrector 
method  and  uses  a  constant  bank  angle  assumption  to  target  the  exit  conditions.  Prior  to 
the  addition  of  the  energy  management  phase,  PredGuid  used  only  the  targeting 
algorithm  throughout  the  trajectory.  Several  heuristic  approximations  were  incorporated 
into  PredGuid  in  order  to  optimize  it  for  a  particular  vehicle  and  mission.  These 
heuristics  and  the  constant  bank  angle  assumption  imposed  severe  limitations  on  the 
usability  of  the  algorithm.  Numerous  improvements  and  changes  were  made  over  the 
course  of  this  study  while  still  preserving  the  basic  intent  of  the  targeting  algorithm.  The 
description  contained  in  this  chapter  describes  the  current  targeting  algorithm  while 
highlighting  major  changes  and  improvements  to  the  original  algorithm. 

Figure  3.4  depicts  the  functional  flow  of  the  targeting  algorithm.  The  following  brief 
outline  summarizes  the  basic  algorithm: 

1 .  Set  the  bank  angle  to  the  current  commanded  bank  angle  {<f>CUD  ) 

2.  Predict  the  final  apogee  using  current  states  and  a  constant  bank  angle  assumption 

3.  Update  the  bank  angle  by  an  amount  A <f>  to  null  the  target  miss 

4.  Repeat  steps  2  and  3  until  guidance  converges  to  the  bank  angle  which  will  hit  the 
target  conditions  within  tolerance,  or  until  the  correction  limit  is  reached 


Figure  3.4:  Targeting  Functional  Flow 
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3.3.1  Predictor-Corrector  Method 

The  predictor-corrector  is  the  main  portion  of  the  targeting  algorithm.  In  the  predictor, 
the  current  states  are  propagated  forward  to  determine  the  resulting  apogee  based  on  a 
constant  bank  angle  assumption.  The  corrector  then  updates  the  bank  angle  guess  in 
order  to  null  the  target  miss.  The  cycle  continues  until  a  bank  angle  is  found  which  meets 
the  apogee  miss  tolerance  or  until  a  maximum  number  of  corrections  is  reached.  The 
number  of  corrections  is  limited  to  improve  running  time  while  still  allowing  guidance  to 
converge  to  a  solution. 


3.3.1. 1  Predictor 


The  numeric  prediction  algorithm  uses  a  Fourth  Order  Runge-Kutta  method  to 
accomplish  the  integration.  The  vector  equations  to  be  integrated  are: 


dr  _  ^ 
dt 


(3.23) 


dv 

—  -  a 
dt 


(3.24) 


The  initial  conditions  are  the  current  position  and  velocity  vectors.  Acceleration  is 
defined  as  the  vector  sum  of  the  gravitational  and  aerodynamic  accelerations: 

a=ag+aa  (3.25) 


The  Earth  gravity  model  used  is  identical  to  the  model  used  in  the  environment  and  is 
described  in  Chapter  2  with  the  final  acceleration  due  to  gravity  described  in  Eq.  2.12  and 
presented  below  as  Eq.  3.26. 
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(3.26) 


Aerodynamic  acceleration  is  further  decomposed  into  acceleration  due  to  drag  and 
acceleration  due  to  lift.  Drag  acts  in  the  opposite  direction  of  the  relative  velocity  and  lift 
acts  perpendicular  to  the  relative  velocity. 

a  a  —  -aJs  ~<*lK  (3-27) 


where 
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_  PVre, 

an  = 

2Cb 

(3.28) 

L 

UL  ~  aD 

(3.29) 

is  =  unit(vrel ) 

(3.30) 

J,  =  unites  x  F) 

(3.31) 

x0cos(^)+Asinfe) 

(3.32) 

<j>p  is  the  bank  angle  used  in  the  prediction.  The  prediction  program  propagates  forward 
in  time  until  the  vehicle  captures  or  escapes  the  atmosphere.  Capture  is  indicated  by  the 
vehicle  dipping  below  200,000  feet  (60.96  km) in  altitude  while  simultaneously  having 
both  negative  altitude  rate  and  acceleration. 

3.3.1. 2  Corrector 

The  correction  algorithm  uses  various  methods  to  compute  the  corrected  bank  angle 
based  on  the  direction  and  severity  of  the  target  misses  from  previous  guesses.  These 
methods  include  interpolation,  extrapolation,  marching  out  of  the  capture  region,  and 
‘smart-guessing’.  The  corrector  functional  flow  is  one  of  the  most  complex  in  the  entire 
program  and  is  depicted  in  Figure  3.5. 
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From  this  figure,  the  six  distinct  methods  of  generating  a  new  bank  angle  guess  are 
visible: 

1 .  Use  previous  bank  angle 

2.  Interpolate  a  high  apogee  bank  angle  and  low  apogee  bank  angle 

3.  Interpolate  between  a  high  apogee  bank  angle  and  a  captured  bank  angle 

4.  Extrapolate  from  two  high  or  two  low  apogee  bank  angles 

5.  March  out  of  the  capture  region 

6.  ‘Smart  guess’  based  on  apogee  miss 

Two  significant  changes  were  made  to  this  portion  of  the  algorithm.  First,  in  the  original 
targeting  algorithm.  Method  6  in  the  list  above  used  a  stored  sensitivity  of  final  apogee  to 
bank  angle  {dRa  /  d<f)\o  determine  the  bank  angle  to  try  based  on  one  ‘good’  guess.  A 
‘good’  guess  is  defined  as  a  guess  not  resulting  in  a  captured  trajectory.  The  stored 
sensitivity  was  specific  to  a  mission  and  vehicle.  Applying  it  to  different  conditions  was 
resulting  in  poor  bank  angle  corrections.  This  method  was  replaced  with  a  simple  ‘smart¬ 
guessing’  algorithm  which  reduced  or  increased  the  bank  angle  by  one  degree  if  the 
predicted  apogee  was  lower  or  higher  than  the  target,  respectively.  While  less 
computationally  efficient  than  use  of  a  stored  sensitivity,  this  method  is  more  generic  and 
was  better  able  to  manage  the  wide  range  of  entry  conditions  examined. 

The  second  change  resulted  from  the  fact  that  the  trajectories  studied  had  a  much  higher 
sensitivity  of  final  apogee  to  bank  angle  change  than  the  original  application  for  which 
the  targeting  algorithm  was  designed.  With  these  high-energy  trajectories,  extremely 
small  bank  angle  changes  yield  large  changes  in  the  predicted  apogee.  Guidance  was 
unable  to  converge  to  a  solution  on  each  execution  because  the  number  of  iterations  was 
insufficient  and  the  minimum  allowable  bank  angle  correction  was  too  small.  The 
number  of  iterations  allowed  per  guidance  cycle  was  increased  and  the  criteria  for  hitting 
the  target  was  loosened  at  high  velocities.  As  the  velocity  decreases,  the  criteria  for 
hitting  the  target  narrows.  Additionally,  the  lower  limit  imposed  on  the  size  of  the  bank 
angle  correction  was  removed  so  that  corrections  could  be  as  small  as  the  sensitivity 
required. 

3.4  Lateral  Control 

Bank  maneuvers  generate  a  component  of  lift  in  the  out-of-plane  direction  causing  out- 
of-plane  position  (dr)  and  out-of-plane  velocity  (<9v)  errors.  Since  bank  angle  is  the  only 
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control  variable,  both  errors  cannot  be  controlled  simultaneously.  For  this  thesis,  the 
desired  orbital  plane  is  assumed  to  be  the  same  as  the  initial  orbital  plane.  The 
perpendicular  to  the  desired  orbital  plane  is  calculated  using  the  initial  position  and 
velocity  vectors  at  El. 

Iyd  =  ,>,,X2"*|  (3-33) 

K*  X  rinit  I 

Out  of  plane  position  and  velocity  are  calculated  by  taking  dot  product  of  the  position  and 
velocity  vectors,  respectively,  with  the  perpendicular. 

(3.34) 

er=r-Iyd  (3.35) 

This  thesis  elects  to  control  out  of  plane  velocity.  If  bank  reversals  are  kept  frequent  and 
even,  position  error  should  be  minimal.  Previously,  in  PredGuid,  a  lateral  corridor  with 
maximum  out  of  plane  velocity  limits  was  determined  by  the  user.  One  change 
incorporated  into  the  targeting  algorithm  was  the  removal  of  heuristic  out  of  plane 
velocity  limits.  They  were  replaced  with  limits  that  are  based  on  a  percentage  of  the 
forward  velocity.  As  a  result,  the  lateral  corridor  still  narrows  with  time,  but  narrows  as  a 
function  of  forward  velocity  rather  than  by  heuristic  limits  which  would  only  work  for 
specific  cases.  Figure  3.6  demonstrates  a  nominal  lateral  corridor. 


Figure  3.6:  Lateral  Corridor  Description 
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3.5  Command  Incorporation 


The  PredGuid  was  designed  to  bias  to  a  neutral  bank  angle  in  order  to  save  command 
margin  for  possible  dispersions  late  in  the  trajectory.  Without  command  biasing, 
dispersions  could  cause  control  saturation  and  subsequent  target  miss.  An  empirical 
velocity  was  chosen  at  which  to  begin  biasing  the  commanded  bank  angle  during 
targeting.  The  bias  feature  remains  in  the  guidance  algorithm,  however  the  velocity 
restriction  was  removed.  The  algorithm  biases  during  the  entire  targeting  phase.  Biasing 
is  unnecessary  during  the  energy  management  phase  and  actually  inhibits  the 
performance  of  that  phase.  In  order  to  avoid  sudden  control  corrections  acceleration 
spikes  at  the  phase  change  point,  the  bias  angle  was  set  to  1 10° .  By  matching  the  phase 
change  criteria,  there  should  be  a  smooth  transition  between  phases. 
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Chapter  4 

Corridor  and  Test  Case  Determination 


4.1  Overview 

As  discussed  in  Chapter  1,  expansion  of  the  capturable  corridor  was  one  of  the  main 
purposes  of  this  thesis.  The  first  step  in  measuring  improvement  is  to  determine  current 
capability.  This  chapter  describes  how  the  theoretical  and  capturable  corridors  are 
determined  and  presents  the  corridors  for  the  lift-to-drag  ratios  and  velocities  examined. 
Furthermore,  this  chapter  describes  initial  conditions  used  in  the  dynamic  simulation. 

4.2  Earth  Corridor  Approximation 

The  theoretical  corridor  is  defined  as  the  range  of  flight  path  angles  through  which  the 
vehicle  is  theoretically  capable  of  accomplishing  successful  aerocapture  where  successful 
aerocapture  is  defined  as  reaching  the  target  apogee  with  no  other  constraints.  The 
shallow  end  of  the  corridor  is  the  flight  path  angle  at  which  the  vehicle  stays  just  below 
the  target  apogee  by  flying  a  full  lift  down  bank  angle  throughout  the  entire  atmosphere. 
The  steep  end  of  the  corridor  is  the  flight  path  angle  at  which  the  vehicle  just  reaches  the 
target  apogee  by  flying  a  full  lift  up  bank  angle  throughout  the  atmosphere.  The  steep 
trajectory  sees  the  greatest  acceleration,  heating  rate,  and  dynamic  pressure  while  having 
a  shorter  duration.  The  shallow  trajectory  has  low  acceleration,  heating  rate,  and 
dynamic  pressure,  but  is  long  in  duration  and  yields  higher  heating  loads.  The  theoretical 
corridor  is  dependent  on  lift-to-drag  ratio.  Figure  4.1  illustrates  the  theoretical  corridor. 

Realistically,  a  vehicle  cannot  fly  the  full  theoretical  corridor  as  defined  above. 
Structural  considerations  limit  the  range  of  flight  path  angles  which  the  vehicle  can 
successfully  fly.  As  already  discussed,  on  the  steep  end  of  the  corridor,  the  vehicle  must 
fly  a  full  lift  up  bank  angle  for  the  entire  trajectory.  If  structural  constraints  limit  the 
maximum  deceleration,  a  more  accurate  description  of  the  steep  end  of  the  corridor 
would  be  the  steepest  flight  path  angle  for  which  the  vehicle  can  fly  full  lift  up  and 
remain  below  an  established  deceleration  limit.  This  thesis  chose  a  maximum  limit  of  10 
g’s,  typical  of  a  robotic  mission. 
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Full  Lift  Up 


Figure  4.1:  Theoretical  Corridor 

The  deceleration  is  calculated  using  Eq.  3.1  and  presented  below  as  Eq.  4.1. 


S © 

Constraints  may  also  be  place  on  heating  rate,  heating  load,  and  AF  to  raise  perigee; 
however,  this  study  did  not  examine  these  constraints.  Flight  path  angles  me'eting  all  the 
constraints  make  up  the  true  theoretical  corridor.  Figure  4.2  through  Figure  4.7 
summarize  the  theoretical  corridors  for  the  lift-to-drag  ratios  examined  as  a  function  of 
hyperbolic  excess  velocity  (Vx).  These  theoretical  corridors  include  only  the 
acceleration  limit  constraints.  The  width  of  the  theoretical  corridor  increases  with 
increased  lift-to-drag  ratio  and  decreasing  Vm . 


44 


Flight  Path  Angle,  y,  (deg) 

Figure  4.2:  Theoretical  Corridor,  L/D  =  0.25 


Figure  4.3:  Theoretical  Corridor,  L/D  =  0.50 


Examination  of  the  acceleration,  dynamic  pressure,  and  heating  rate  over  time 
demonstrates  the  large  difference  between  trajectories  on  the  opposite  ends  of  the 
corridor.  Figure  4.8  through  Figure  4.10  demonstrate  this  concept  for  a  representative 
sample  case  with  lift-to-drag  ratio  of  0.75  and  Vx  =  6  km/ s  .  For  this  case,  the  heating 
load  on  the  shallow  end  of  the  corridor  was  24,214  Btu/ft2  (275,199  kJ/m2)  and  the 
heating  load  on  the  steep  end  of  the  corridor  was  18,062  Btu/ft2  (204,984  kJ/m2). 


Time  (s) 

Figure  4.8:  Acceleration  vs.  Time,  Corridor  Extremes 
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4.3  Test  Case  Determination 

Entry  flight  path  angles  were  determined  using  basic  astrodynamics.  Six  values  are 
required  to  fully  describe  an  orbit  and  the  vehicle  position  on  the  orbit.  The  simulation 
was  initialized  with  position  and  velocity  vectors  derived  from  these  values.  The  six 
values  chosen  for  this  thesis  were  hyperbolic  excess  velocity  (Vx),  altitude  at  El  (a£/), 
inclination  (/),  vacuum  periapse  altitude  (hp),  argument  of  perigee  (co),  and  longitude  of 
the  ascending  node(fi).  Four  of  these  values  remained  constant  for  all  cases: 

aEi  ~  400,000  ft  (121.92  km) 

0)  =  0° 

Q  =  0° 
i  =  5° 

Altitude  at  El  was  chosen  based  on  the  top  of  the  sensible  atmosphere,  standardized  for 
Earth  at  400,000  ft.  Argument  of  perigee  and  longitude  of  the  ascending  node  were 
chosen  for  simplicity.  Inclination  was  chosen  to  ensure  that  argument  of  perigee  and 
longitude  of  the  ascending  node  would  be  defined. 

Reference  [3]  compares  hyperbolic  excess  velocities  for  other  planets  and  moons  and 
translates  them  to  equivalent  Earth  demonstration  velocities.  Earth  approach  velocities 
which  effectively  simulate  approaches  to  other  planets  and  moons  range  from  9  to  14 
km/s  (29528  to  45931  ft/s).  The  velocities  examined  in  this  thesis  were  1,  3,  6,  10,  14, 
and  18  km/s  (3281,  9843,  19685,  32808,  45931,  and  59055  ft/s)  which  completely 
encompasses  the  Earth  demonstration  velocity  range  from  reference  [3]  and  extends 
beyond  it  for  further  examination.  At  each  velocity,  vacuum  periapse  altitudes  were 
chosen  to  yield  the  desired  flight  path  angles.  The  initial  conditions  were  independent  of 
lift-to-drag  ratio. 
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Chapter  5 
Results 


5.1  Overview 

This  chapter  presents  the  results  of  the  algorithm  enhancements.  The  first  section 
provides  verification  that  the  energy  management  guidance  is  working  correctly. 
Capturable  corridor  results  are  presented  in  the  next  section.  Finally,  the  third  section 
discusses  various  metrics  investigated  in  an  effort  to  analyze  the  effectiveness  of  the 
enhanced  algorithm. 

5.2  Guidance  Verification 

As  discussed  in  Chapter  3,  the  energy  management  phase  seeks  to  maintain  a  constant 
cruise  condition  until  sufficient  energy  has  been  depleted  to  allow  successful  targeting. 
Select  cases  were  examined  to  ensure  correct  performance.  For  a  nominal  case,  several 
observations  would  be  expected.  The  vehicle  would  attempt  to  drive  the  altitude  rate  to 
zero  by  initially  flying  full  lift  up.  After  the  altitude  rate  settled  to  zero,  the  bank  angle 
would  need  to  be  lift  down  to  keep  the  vehicle  in  the  atmosphere.  However,  the  vehicle 
will  slowly  push  more  lift  up  to  generate  enough  lift  to  maintain  altitude  as  velocity 
decreases.  Acceleration  would  initially  peak  as  the  vehicle  attempts  to  null  the  altitude 
rate.  Finally,  drag  attempts  to  follow  the  reference  drag  curve.  Figure  5.1  through  Figure 
5.5  depict  a  nominal  case  ( Vx=\4km/s ,  L/D-0.15,  y  =  5.6°)  demonstrating  the 
energy  management  phase  behaves  in  this  manner  prior  to  phase  change.  Note  that  due 
to  the  decoupling  of  vertical  and  lateral  guidance  and  the  assumption  of  instantaneous 
bank  angle  changes,  bank  angle  results  are  presented  as  the  absolute  value  of  the 
commanded  bank  angle,  \</>c  \ .  This  method  allows  easier  observation  of  bank  angle 
values  and  trends  without  losing  vertical  guidance  information. 
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Figure  5.4:  Commanded  Bank  Angle 
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Figure  5.5:  Drag  Acceleration  vs.  Time,  Voo  =  14  km/s,  L/D  =  0.75,  y  =  6.9° 

Figure  5.6  through  Figure  5.19  depict  commanded  bank  angle  histories  for  various 
velocities  and  lift-to-drag  ratios.  Each  figure  depicts  a  velocity  and  lift-to-drag  ratio  pair. 
The  lift-to-drag  ratios  depicted  are  for  the  low  (0.25),  middle  (0.75),  and  high  (1.5)  values 
in  the  range  studied.  All  velocities  for  which  guidance  was  able  to  capture  a  portion  of 
the  corridor  are  graphed  for  each  lift-to-drag  ratio.  Each  figure  plots  bank  angle 
trajectories  for  the  steepest  and  shallowest  flight  path  angles  captured  as  well  as  a  flight 
path  angle  in  the  middle  between  these  extremes. 
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Figure  5.8:  Commanded  Bank  Angle  vs.  Time,  Voo  =  6  km/s,  L/D  =  0.25 
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Figure  5.9:  Commanded  Bank  Angle  vs.  Time,  V*  =  10  km/s,  L/D  -  0.25 
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Figure  5.10:  Commanded  Bank  Angle  vs.  Time,  V*,  =  1  km/s,  L/D  =  0.75 
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Figure  5.11:  Commanded  Bank  Angle  vs.  Time,  V*  =  3  km/s,  L/D  =  0.75 
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Figure  5.12:  Commanded  Bank  Angle  vs.  Time,  V*  =  6  km/s,  L/D  =  0.75 
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Figure  5.13:  Commanded  Bank  Angle  vs.  Time,  Voo  =  10  km/s,  L/D  -  0.75 
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Figure  5.14:  Commanded  Bank  Angle  vs.  Time,  Voo  =  14  km/s,  L/D  =  0.75 
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Figure  5.15:  Commanded  Bank  Angle  vs.  Time,  V«,  =  1  km/s,  L/D  =  1.5 
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Figure  5.16:  Commanded  Bank  Angle  vs.  Time,  V«>  =  3  km/s,  L/D  -  1.5 
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Figure  5.17:  Commanded  Bank  Angle  vs.  Time,  V®  =  6  km/s,  L/D  =  1.5 
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Figure  5.18:  Commanded  Bank  Angle  vs.  Time,  V»  =  10  km/s,  L/D  =  1.5 
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Figure  5.19:  Commanded  Bank  Angle  vs.  Time,  Voo  =  14  km/s,  L/D  =  1.5 
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These  figures  show  several  important  trends.  First,  at  low  velocities,  for  all  lift-to-drag 
ratios,  the  guidance  transitions  to  targeting  almost  immediately.  This  trend  demonstrates 
that  the  phase  check  is  working  properly.  Lower  velocities  do  not  require  as  much 
energy  depletion  for  successful  targeting  and  should  transition  to  targeting  earlier. 

Second,  at  higher  velocities,  shallow  flight  path  angles  require  longer  energy 
management  phases  prior  to  transition  to  targeting.  As  lift-to-drag  ratio  increases  at  high 
velocities,  regardless  of  flight  path  angle,  the  vehicle  also  spends  longer  in  the  energy 
management  phase.  These  trends  demonstrate  that  the  energy  management  phase  and 
phase  change  logic  are  working  correctly.  At  shallow  flight  path  angles,  the  vehicle  does 
not  dig  as  deeply  into  the  atmosphere.  Vehicles  with  high  lift-to-drag  ratios  null  the 
altitude  rate  more  quickly  so  they  also  do  not  dig  as  deeply  into  the  atmosphere.  Due  to 
lower  drag  at  these  atmospheric  levels,  the  vehicle  must  spend  longer  in  the  energy 
management  phase  to  deplete  the  same  A  V . 

Third,  almost  all  of  the  trajectories  have  a  ripple  in  the  commanded  bank  angle  histories 
during  later  portions  of  the  targeting  phase.  This  ripple  corresponds  to  the  apogee  miss 
criteria  changeover  point.  At  this  point  in  targeting,  the  apogee  miss  criteria  for  a  ‘good’ 
solution  is  reduced  from  40  nm  (74.08  km)  to  3  nm  (5.57  km).  This  trend  demonstrates 
that  the  targeting  guidance  is  successfully  converging  to  bank  angle  solutions.  When  the 
acceptance  criteria  changes,  guidance  immediately  takes  action  to  hone  in  on  the  target. 

5.3  Capturable  Corridor  Coverage 

As  discussed  in  Chapter  4,  the  capturable  corridor  is  defined  as  the  portion  of  the 
theoretical  corridor  for  which  the  algorithm  is  able  to  successfully  reach  the  target 
without  exceeding  a  maximum  acceleration  of  10  g’s.  Figure  5.20  through  Figure  5.25 
show  how  the  enhanced  and  original  algorithms’  capturable  corridors  compare  to  the 
theoretical  corridors  presented  in  Chapter  4.  Corridors  are  depicted  for  lift-to-drag  ratios 
as  a  function  of  Vm . 
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Flight  Path  Angle,  y,  (deg) 

Figure  5.20:  Theoretical  and  Capturable  Corridors,  L/D  =  0.25 


Flight  Path  Angle,  y,  (deg) 


Figure  5.21:  Theoretical  and  Capturable  Corridors,  L/D  =  0.50 
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Flight  Path  Angle,  y,  (deg) 

Figure  5.24:  Theoretical  and  Capturable  Corridors,  L/D  =  1.25 


Flight  Path  Angle,  yf  (deg) 

Figure  5.25:  Theoretical  and  Capturable  Corridors,  L/D  =  1.50 


As  you  can  see  from  the  plots,  the  algorithm  enhancements  greatly  improved  corridor 
coverage,  especially  for  the  14  km/s  cases.  Figure  5.26  shows  the  percent  improvement 
over  PredGuid  of  corridor  captured  for  the  various  lift-to-drag  ratios  and  values  of  Vx . 


Lift-to-Drag  Ratio 

Figure  5.26:  Percent  Improvement  of  Corridor  Captured 

Most  of  the  corridor  coverage  gain  occurred  on  the  steep  end  of  the  corridor.  Also,  two 
distinct  areas  of  the  theoretical  corridor  lacked  coverage:  the  shallow  end  of  the  corridor 
(much  of  which  lost  area  that  PredGuid  was  able  to  cover)  and  the  steep  flight  path  angles 
with  low  velocity.  These  results  can  be  explained  by  the  nature  of  how  each  algorithm 
works. 

5.3.1  Steep  Flight  Path  Angle  Corridor  Gain 

Steep  end  corridor  coverage  gain  is  the  result  of  the  enhanced  algorithm’s  ability  to 
remain  at  a  particular  altitude  as  long  as  necessary  to  remove  energy.  The  original 
algorithm  could  only  penetrate  deeper  into  the  atmosphere  to  generate  more  drag  and 
deeper  vehicle  atmospheric  penetration  yields  high  acceleration  forces.  Figure  5.27 
through  Figure  5.29  depict  an  example  case  (F*,  =  6  Atm/s  ,  L/D  =  0.75  ,  and  y  =  5.9°) 
for  which  both  algorithms  are  able  to  reach  the  target  apogee,  but  the  original  algorithm 
causes  a  much  higher  (and  unacceptable)  acceleration. 
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Figure  5.29:  Commanded  Bank  Angle  vs.  Time,  V,*,  =  6  km/s,  L/D  =  0.75,  y  -  5.9° 

5.3.2  Enhanced  Algorithm  Corridor  Deficiencies 

There  are  two  distinct  areas  of  deficiency:  the  shallow  end  of  the  corridor  and  the  steep 
flight  path  angles  with  low  hyperbolic  excess  velocity. 

5.3.3  Shallow  End  Corridor  Loss 

On  the  shallow  end  of  the  corridor,  the  vehicle  needs  to  fly  more  lift  down  throughout  the 
trajectory  in  order  to  reach  the  target.  However,  the  enhanced  algorithm  automatically 
attempts  to  first  null  the  altitude  rate  by  flying  full  lift  up  before  checking  for  a  phase 
change.  By  the  time  the  vehicle  flips  to  full  lift  down,  the  altitude  rate  cannot  be  slowed 
and  the  vehicle  skips  out  of  the  atmosphere.  Alternatively,  PredGuid  allowed  the  vehicle 
to  fly  lift  down  for  the  entire  trajectory.  As  a  result,  the  original  algorithm  had  better 
shallow  flight  path  angle  coverage.  Figure  5.30  through  Figure  5.32  depict  the  bank 
angle,  altitude  and  altitude  rate  profiles  for  an  example  case  (Vx  =  3  km/s ,  L/D  -  0.75 , 
and  y  =  4.75° )  in  which  PredGuid  was  able  to  target  the  final  apogee  but  the  enhanced 
algorithm  was  not.  The  enhanced  algorithm  never  transitioned  out  of  the  energy 
management  phase  and  was  unable  to  remain  in  the  atmosphere  even  flying  full  lift  down 
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because  the  vehicle  did  not  penetrate  far  enough  into  the  atmosphere  to  generate  enough 
lift  to  null  the  altitude  rate. 
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Figure  5.32:  Altitude  Rate  vs.  Time,  V,*,  =  3  km/s,  L/D  —  0.75,  y  —  4.75' 


Figure  5.33:  Acceleration  vs.  Time,  Vo,  =  3  km/s,  L/D  =  0.75,  y  -  4.75° 


This  same  trend  continued  throughout  the  higher  hyperbolic  excess  velocities.  Even  the 
vehicles  with  the  highest  lift-to-drag  ratios  were  still  unable  to  capture  on  the  shallow 
side  of  the  corridor.  The  ‘shallow’  range  extended  to  steeper  and  steeper  flight  path 
angles  as  Vx  increased.  Figure  5.34  and  Figure  5.35  demonstrate  a  sample  case  with  a 
fairly  steep  flight  path  angle  but  high  velocity  and  low  lift-to-drag  ratio  (Vn  =14  km/ s , 
L/D  =  0.25,  y  =  6.5°).  Figure  5.36  and  Figure  5.37  show  a  sample  case  with  high 
velocity  and  high  lift-to-drag  ratio  (Vx  =18  km/ s  ,  LjD  =  1.5,  y  =  6.25°  ). 


Time  (s) 

Figure  5.34:  Commanded  Bank  Angle  vs.  Time,  V*,  =  14  km/s,  L/D  =  0.25,  y  =  6.5° 
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Figure  5.37:  Altitude  Rate  vs.  Time,  Vo,  =  18  km/s,  L/D  =  1.5,  y  =  6.25° 

As  you  can  see,  the  vehicle  response  is  virtually  identical  for  all  of  the  cases  shown.  For 
shallow  entries,  flying  full  lift  up  set  up  orbital  conditions  from  which  the  vehicle  could 
not  recover.  The  vehicle  never  transitioned  out  of  the  energy  management  phase  and  was 
unable  to  remain  in  the  atmosphere  despite  flying  full  lift  down. 

5.3.4  Steep  Flight  Path  Angles  with  Low  Velocity 

Figure  5.38  through  Figure  5.40  show  a  sample  case  with  a  steep  flight  path  angle  beyond 
the  capturable  corridor  of  the  enhanced  algorithm  ( Vm  =  6  km  / s ,  L/D  =  0.15, 
y  =  7.25°  ). 

The  algorithm  successfully  targeted  the  apogee,  but  the  acceleration  loads  were 
unacceptably  high.  The  peak  acceleration  occurs  after  the  vehicle  has  already 
transitioned  to  targeting.  At  these  steep  flight  path  angles,  transition  to  targeting  occurred 
too  early  before  guidance  was  able  to  null  the  altitude  rate  and  dissipate  energy. 
Transition  occurred  early  because  the  predicted  trajectories  indicated  capture.  This 
problem  only  occurs  at  low  velocities  which  are  more  susceptible  to  capture.  Once 
guidance  switched  to  the  targeting  phase  prematurely,  the  excessive  acceleration 
problems  that  plagued  PredGuid  resurfaced. 
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Figure  5.40:  Altitude  Rate  vs.  Time,  Vx  =  6  km/s,  L/D  =  0.75,  y  =  7.25° 

5.4  Algorithm  Metrics 

Improved  corridor  coverage  should  indicate  a  better  algorithm.  However,  a  discussion  of 
improvement  based  solely  on  corridor  coverage  is  incomplete  for  this  thesis  because  only 
acceleration  loads  were  considered  in  the  corridor  constraints.  There  are  several  other 
factors  to  consider  when  evaluating  performance  including  heating  rate,  heating  load,  and 
AV  to  raise  perigee.  Each  of  these  metrics  as  well  as  acceleration  will  be  discussed  in 
this  section. 

5.4.1  Acceleration 

Figure  5.41  shows  maximum  vehicle  acceleration  as  a  function  of  flight  path  angle  for 
the  original  and  enhanced  algorithms,  respectively.  Each  dot  represents  a  lift-to-drag 
ratio,  Vm,  and  flight  path  angle  trio.  As  you  can  see  from  the  graphs,  the  original 
algorithm  yielded  much  higher  maximum  accelerations  than  the  enhanced  algorithm  at 
the  same  flight  path  angles.  These  results  are  not  surprising  because  the  algorithm  was 
designed  to  operate  at  a  lower  acceleration  over  a  longer  period  of  time.  These  figures 
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also  help  visualize  the  expanded  corridor  coverage.  If  the  acceleration  limit  were  placed 
at  a  different  level  than  10  g’s,  the  capturable  corridor  limits  would  change. 

PredGuid  Enhanced 


Figure  5.41:  Maximum  Acceleration  vs.  Flight  Path  Angle 


5.4.2  Heating  Rate 


Maximum  heating  rate  is  another  factor  that  must  be  considered.  Again,  at  similar  flight 
path  angles,  the  enhanced  algorithm  had  lower  maximum  heating  rates  as  shown  in 
Figure  5.42.  The  lower  maximum  heating  rates  are  a  product  of  the  lower  decelerations 
and  shallower  trajectories.  The  few  high  heating  rates  visible  in  the  enhanced  algorithm 
plot  are  high  velocity  trajectories.  There  are  no  figures  to  compare  them  against  because 
the  original  algorithm  does  not  work  for  these  cases. 

PredGuid  Enhanced 


y  (deg) 


Y  (deg) 


Figure  5.42:  Maximum  Heating  Rate  vs.  Flight  Path  Angle 
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5.4.3  Heating  Load 


The  enhanced  algorithm  typically  resulted  in  longer  periods  of  time  spent  in  the 
atmosphere  which  would  normally  cause  higher  heating  loads.  However,  as  discussed 
above,  the  heating  rates  were  lowered.  The  longer  atmospheric  flights  do  not  appear  to 
cause  higher  heating  loads  due  to  the  decreased  heating  rates.  Heating  loads  as  a  function 
of  time  are  presented  in  Figure  5.43.  The  few  high  heating  loads  visible  in  the  enhanced 
algorithm  graph  are  for  high  velocity  trajectories.  As  before,  these  cannot  be  compared 
to  known  data  because  the  original  algorithm  failed  to  capture  any  portion  of  the  corridor 
at  those  velocities. 


Figure  5.43:  Heating  Load  vs.  Flight  Path  Angle 


5.4.4  Velocity  Change  to  Raise  Perigee 

The  final  factor  examined  was  the  AV  required  to  circularize  the  orbit  after  aerocapture. 
As  seen  in  Figure  5.44,  AV  required  to  circularize  held  fairly  stable  for  all  flight  path 
angles  because  both  PredGuid  and  the  enhanced  algorithm  used  the  same  targeting 
algorithm.  Even  the  flight  path  angles  on  the  portion  of  the  corridor  captured  by  the 
enhanced  algorithm  still  had  approximately  the  same  AV  as  the  portion  captured  by  the 
original  algorithm.  PredGuid  is  an  excellent  targeting  algorithm  and  the  energy 
management  additions  did  not  inhibit  its  function. 
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Figure  5.44:  AV  to  Circularize  vs.  Flight  Path  Angle 


5.5  Lateral  Control 

The  updated  lateral  control  portion  of  the  guidance  program  worked  fairly  well.  Figure 
5.45  depicts  the  lateral  corridor  information  for  a  nominal  case  which  remained  in  the 
energy  management  phase  for  a  significant  period  of  time  (Vx  =14  km  I  s ,  L/D  =  0.75  , 
7  =  6.9°). 


Figure  5.45:  Out  of  Plane  Velocity  vs.  Time,  V*,  =  14  km/s,  L/D  =  0.75,  y  -  6.9' 


As  you  can  see  from  the  above  figure,  the  bank  angle  reversals  are  executed  in  a  timely 
manner.  The  use  of  percentage  of  forward  velocity  allowed  a  gradual  narrowing  of  the 
corridor  early  in  the  trajectory  and  prevented  excessive  reversals.  However,  the  corridor 
did  not  narrow  as  quickly  near  the  end  of  the  trajectory  because  the  vehicle  is 
decelerating  at  a  slower  rate.  For  higher  lift-to-drag  ratio  cases,  bank  reversals  occurred 
more  frequently  because  higher  lift-to-drag  ratios  cause  more  out  of  plane  lift.  Figure 
5.46  depicts  corridor  information  for  a  higher  lift-to-drag  ratio  case.  Finally,  a  case  was 
examined  which  did  not  remain  in  the  energy  management  phase  for  an  extended  period 
(Vm  =  6  km/s ,  L/D  =  0.25 ,  y  =  6.0° ).  Figure  5.47  depicts  these  results.  Even  without 
the  gradual  narrowing  of  the  lateral  corridor  during  energy  management,  the  guidance 
still  managed  to  keep  the  lateral  error  within  the  corridor. 


Figure  5.46:  Out  of  Plane  Velocity  vs.  Time,  V*  =  14  km/s,  L/D  =  1.5,  y  =  6.9° 
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Figure  5.47:  Out  of  Plane  Velocity  vs. 
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Chapter  6 
Conclusions 

6.1  Summary  and  Conclusions 

This  thesis  sought  to  enhance  an  existing  numerical  predictor-corrector  aerocapture 
guidance  algorithm.  The  enhancements  were  accomplished  by  implementing  an  energy 
management  phase  prior  to  targeting  exit  conditions  and  replacing  heuristic  features  with 
more  generic  features.  The  purpose  of  the  energy  management  phase  is  to  deplete 
sufficient  energy  to  allow  targeting  of  the  exit  conditions  in  the  targeting  phase.  The 
vehicle  response  during  energy  management  can  be  modeled  as  a  second-order 
spring/mass/damper  system.  The  bank  angle  is  determined  by  calculating  the  bank  angle 
for  zero  altitude  acceleration  and  then  adding  an  altitude  rate  damper  and  a  term  to  null 
the  drag  profile  deviation.  Phase  change  occurs  when  two  conditions  have  been  met: 
First,  if  at  the  given  point  in  the  trajectory,  the  vehicle  could  fly  a  constant  bank  angle  of 
1 10°  for  the  remainder  of  the  trajectory  and  have  the  resulting  apogee  be  below  or  within 
a  given  tolerance  above  the  target  apogee.  Second,  the  predicted  final  energy  must 
indicate  that  the  vehicle  would  be  on  an  elliptical,  not  hyperbolic,  trajectory. 

Implementation  of  the  energy  management  phase  allowed  removal  of  several  heuristic 
features.  Modeling  of  a  separate  lift  down  phase  was  replaced  by  using  a  lift  down 
condition  to  determine  phase  change  and  biasing  to  the  same  lift  down  condition  during 
targeting.  Use  of  a  heuristic  sensitivity  to  calculate  the  first  corrected  bank  angle  was 
replaced  by  a  simple  ‘smart  guessing’  algorithm  based  on  whether  the  predicted  apogee 
was  above  or  below  the  target.  Heuristic  lateral  corridor  boundaries  were  replaced  by 
boundaries  based  on  percentage  of  forward  velocity. 

An  analysis  of  the  resulting  flight  path  angle  entry  corridor  revealed  that  the  enhanced 
algorithm  allowed  1 14%  corridor  coverage  improvement  the  lowest  lift-to-drag  ratio,  and 
379%  for  the  highest  lift-to-drag  ratio.  The  mean  improvement  was  261%  over  the 
PredGuid  capturable  corridor  for  all  velocities  and  lift-to-drag  ratios.  Most  of  the  gain 
occurred  for  higher  flight  path  angles.  However,  results  also  indicated  that  at  the  same 
conditions,  the  enhanced  algorithm  yielded  lower  maximum  accelerations.  Comparable 
heating  rates,  heating  loads,  and  AV  to  raise  perigee  were  observed.  The  AV  to  raise 
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perigee  results  strongly  indicate  that  PredGuid  is  an  excellent  targeting  algorithm  and  the 
enhancements  did  not  inhibit  the  algorithm’s  targeting  performance. 

6.2  Recommendations  for  Future  Work 

There  are  several  areas  of  this  project  which  could  be  examined  in  the  future. 

Shallow  Flight  Path  Angles:  One  of  the  most  important  areas  for  future  study  would  be 
to  examine  how  to  make  the  enhanced  algorithm  work  for  the  shallow  end  of  the  entry 
corridor.  The  original  guidance  algorithm  was  able  to  capture  within  this  region  but  the 
enhanced  algorithm  could  not.  There  should  be  a  recourse  which  enables  the  algorithm 
to  vary  the  initial  guidance  command  rather  than  simply  flying  full  lift  up  to  null  the 
altitude  rate. 


Steep  Flight  Path  Angles  with  Low  Velocity:  This  area  was  the  only  other  ‘hole’  in  the 
enhanced  algorithm’s  corridor  coverage.  Guidance  transitions  between  the  energy 
management  and  targeting  phases  prematurely  at  these  conditions.  Another  condition 
needs  to  be  implemented  into  the  phase  check  to  prevent  this  premature  transition. 


Apogee  Miss  Criteria:  The  only  true  heuristic  remaining  in  the  program  is  the  apogee 
error  miss.  Due  to  the  extreme  sensitivity  of  final  apogee  to  bank  angle  at  high  velocities, 
placing  a  constraint  to  narrow  the  apogee  miss  as  velocity  decreased  was  unavoidable. 
This  constraint  caused  reactionary  spikes  in  almost  every  bank  angle  profile  as  the 
apogee  miss  suddenly  dropped  from  40  nm  to  3  nm  at  a  velocity  of  26,600  ft/s.  A  more 
generic  method,  possibly  based  on  percentage  of  velocity  or  energy  decrease,  to  generate 
a  more  gradual  narrowing  of  the  apogee  miss  should  be  investigated. 


Modeling  Bank  Changes:  This  thesis  assumed  instantaneous  bank  angle  changes.  For 
small  changes  in  bank  angle,  this  assumption  is  fairly  accurate.  However,  as  bank 
changes  increase,  the  reversals  and  corrections  can  take  a  significant  amount  of  time. 
Incorporation  of  actual  vehicle  roll  rate  would  be  more  realistic.  Modeling  should 
include  roll  rate,  roll  acceleration  and  deceleration,  and  the  ability  to  roll  over  or  under 
through  the  shortest  distance. 
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Lateral  Corridor:  The  lateral  control  presented  in  this  thesis  works  reasonably  well. 
However,  there  is  still  significant  out  of  plane  velocity  error  present  at  the  end  of  the 
trajectory.  Another  method  should  be  investigated  which  allows  additional  narrowing  of 
the  lateral  corridor  towards  the  end  of  the  trajectory. 


Robustness  analysis:  This  thesis  assumed  perfect  navigation  and  control.  A  robustness 
analysis  including  density  dispersions,  lift-to-drag  ratio  estimate  errors,  control  errors, 
etc.  would  provide  a  better  metric  of  how  well  the  algorithm  works. 
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Appendix  A 
PredGuid 

A.l  Aerocapture  Executive 

Aerocapture  guidance  uses  a  predictor-corrector  method  to  determine  the  periodic  bank 
angle  commands  that  control  the  spacecraft  to  a  target  apogee.  Bank  angle  commands 
are  given  to  produce  desired  exit  conditions.  Bank  angle  commands  generate  an 
excessive  out  of  plane  velocity;  as  a  result,  periodic  roll  reversals  are  commanded  to  keep 
the  spacecraft  within  the  desired  orbital  plane. 


The  logical  flow  of  the  Aerocapture  Executive  algorithm  is  illustrated  in  Figure  A.l. 
Functional  flow  of  the  Aerocapture  Executive  is  illustrated  in  Figure  A.2 
The  data  inputs  to  pred_guid.m  are  described  in  Table  A-l . 

The  data  outputs  of  pred_guid.m  are  described  in  Table  A-2. 

The  constants  used  in  pred  guid.m  are  described  in  Table  A-3. 


1.  If  this  is  the  first  guidance  pass  (guidance  initialization  flag  is  set  to  one),  execute 
steps  a  through  c.  Otherwise,  proceed  directly  to  step  2: 


a.  Define  the  constants  structure,  c. 


b.  Initialize  the  local  guidance  parameters  structure,  1. 

c.  Execute  the  Guidance  Initialization  procedure  to  fill  1. 

2.  If  the  spacecraft  is  in  the  sensible  atmosphere  (indicated  by  a  user-defined 
constant  g  value),  update  coefficient  of  drag,  density  gain,  and  lift-to-drag  ratio  by 
executing  the  Aerodynamic  Properties  procedure.  Otherwise,  proceed  directly  to 
step  3. 
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3.  If  the  spacecraft  is  in  the  sensible  atmosphere  (indicated  by  a  constant  g  value 
higher  than  the  constant  g  value  required  to  execute  the  Aerodynamic  Properties 
procedure  ensuring  that  the  Aerodynamic  Properties  procedure  filter  has  had 
sufficient  time  to  provide  good  aerodynamic  data),  or  if  this  is  the  first  pass 
(guidance  initialization  flag  is  set  to  one),  then  execute  steps  a  and  b.  Otherwise, 
proceed  directly  to  step  4. 

a.  Set  the  guidance  initialization  flag  to  zero  to  ensure  that  the  variables  are 
not  initialized  again: 

F.  =  0 

ig 

b.  If  the  guidance  pass  number  is  zero,  then  get  desired  bank  angle  from  the 
Predictor-Corrector  Sequencer  procedure.  A  guidance  pass  number  equal 
to  zero  indicates  either  the  first  pass  or  that  you  have  reached  the  limit  of 
guidance  passes  without  updating  the  bank  angle.  The  aerocapture 
guidance  algorithm  is  intended  to  be  run  at  a  frequency  of  1  Hz  to  give  a 
bank  angle  command  once  per  second;  however,  bank  angle  updates  need 
not  be  made  every  second.  The  guidance  pass  limit  variable  enables  the 
user  to  determine  the  frequency  at  which  to  update  the  bank  angle. 

4.  Execute  the  Lateral  Guidance  procedure  to  get  the  sign  of  the  bank  angle. 

5.  Execute  the  Bank  Angle  Command  Incorporation  procedure  to  get  the 
commanded  bank  angle. 

6.  Increment  the  guidance  pass  number: 

GP  =  GP  + 1  (A.1) 

1.  If  you  have  reached  the  guidance  pass  limit  (a  user-defined  number  which 
determines  the  frequency  at  which  to  update  the  bank  angle),  reset  the  guidance 
pass  number  to  zero: 

GP  =  0 
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* 

Command 

Incorporation 


K Commanded 

No - ►/  Bank  Angle  to 

V  Guidance 


Figure  A.l:  Aerocapture  Guidance  Logical  Flow:  pred_guid.m 


Figure  A.2:  Aerocapture  Guidance  Functional  Flow:  pred_guid.m 
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Table  A-l:  Inputs,  predguid.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

i 

A_NAV 

anav 

Vehicle  inertial 
aerodynamic 
acceleration 
vector 

Ft/s2 

double[3] 

Navigation 

i 

ALT_NAV 

^ nav 

Vehicle  altitude 

Ft 

double 

Derived 

i 

G  LOAD 

g 

Vehicle 
aerodynamic 
load  factor 

double 

Derived 

i 

IYD 

I  YD 

Unit  normal  to 
the  desired  orbit 
plane 

double[3] 

Derived 

i 

QJBTOI 

Qbi 

Vehicle  body- 
to-inertial 
attitude 
quaternion 

(-) 

double  [4] 

Navigation 

i 

R_NAV 

Rna, 

Vehicle  inertial 
position  vector 

ft 

double[3] 

Navigation 

i 

V_NAV 

vnav 

Vehicle  inertial 
velocity  vector 

ft/s 

double[3] 

Navigation 

i 

V_NAV_MAG 

\Ka\ 

Vehicle  inertial 
velocity  vector 
magnitude 

ft/s 

double 

Derived 

i 

VRELMAG 

fre,  | 

Vehicle  relative 
velocity  vector 
magnitude 

ft/s 

double 

Derived 

i 

V_REL_NAV 

Ke, 

Vehicle  relative 
velocity  vector 

ft/s 

double[3] 

Derived 

Table  A-2:  Outputs,  pred  guid.m 


Structure 

Name 

Variable 

Name 

Symbol 

Data  Type/ 
Precision 

Destination 

PHICMD 

<t>cmd 

Commanded 
bank  angle 

deg 

double 

Guidance 

System 
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Table  A-3:  Constants,  pred_guid.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

c 

AEROPROPERTYGLOAD 

&  aero 

Minimum  g 
load 

required  to 
update 
aerodynamic 
properties 

G’s 

double 

0.05 

c 

GRUNGUIDANCE 

&  guid 

Minimum  g 
load 

required  to 
run 

guidance 

G’s 

double 

0.075 

c 

GUIDPASSLIM 

r 

° gpl 

Max  number 
of  guidance 
passes 
allowed 
without 
updating 
bank  angle 

(-) 

double 

10 

A.2  Guidance  Initialization 

The  initial_guid.m  function  initializes  the  local  guidance  parameters  structure,  1,  on  the 
first  guidance  pass. 


The  logical  flow  of  the  Guidance  Initialization  algorithm  is  illustrated  in  Figure  A.3. 

The  data  inputs  to  initialguid.m  are  described  in  Table  A-4. 

The  local  guidance  parameters  used  or  updated  in  initial_guid.m  are  described  in  Table 
A-5. 

The  constants  used  in  initial  guid.m  are  described  in  Table  A-6. 


1 .  Set  the  maximum  and  minimum  cosine  of  bank  angle  values  (corresponding  to  the 
minimum  and  maximum  bank  angles,  respectively)  equal  to  the  user-defined 
constant  values 
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C0SWmin  cos(^max  C  r  ) 
cos(^Iiax  =  cos(^min  Cd->r  ) 

2.  Initialize  the  bank  angle  sign  parameters  by  executing  steps  a  through  c: 


a.  Calculate  the  sign  of  the  bank  angle  using  the  currently  velocity  vector  and 
unit  normal  to  the  desired  plane 

sign{<f>)  =  sign(Vnav  x  IYD )  (A.2) 


b.  Initialize  the  commanded  bank  angle  using  this  sign  of  bank  and  the  user- 
defined  initial  value 

=sign{<f\<pEi\  (A-3) 


c.  Initialize  the  previous  sign  of  bank  value  to  the  current  sign  of  bank  value 
siSnWpasl  =sign{(f) 

3.  Initialize  the  cosine  of  the  desired  bank  angle  to  the  user-defined  initial  value 

C0S(^L  =  cos(^L„„ 

4.  Initialize  the  roll  reversal  parameters  by  executing  steps  a  through  c : 


a.  Calculate  the  corridor  slope.  Essentially  we  are  creating  an  equation  of  a  line 
where  the  corridor  slope  and  the  minimum  corridor  out  of  plane  velocity 
define  the  slope  and  the  y  intercept,  respectively.  This  is  a  user-defined 
‘acceptability’  criterion  to  determine  if,  based  on  the  out  of  plane  velocity,  the 
spacecraft  is  within  the  desired  orbital  plane. 


S 


corr 


CL 


c 


cv 


■cv 


(A.4) 


b.  Initialize  the  lift  down  reversal  flag  to  one  (Feature  not  currently  used) 

F  =  1 

r  LDR  1 

c.  Initialize  the  over/under  bank  angle  to  the  predicted  angle  (Feature  not 
currently  used) 

fiou  =  tftou  pred 

5.  Initialize  the  lift  down  phase  variables  by  executing  steps  a  and  b: 

a.  Initialize  the  lift  down  modeling  flag.  A  flag  value  of  zero  deactivates  lift 
down  modeling  while  a  value  of  one  activates  it. 

FMLd  =  Fuldm, 
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b.  Set  the  initial  sine  and  cosine  values  of  the  lift  down  bank  angle  to  the  user- 
defined  initial  values 

COS {</>LD  )  =  COS Cp-tR  ) 

sin(^D)  =  sin(^DmCD^) 

6.  Initialize  aerodynamic  parameters  to  user-defined  initial  values: 


a.  Drag  coefficient 

b.  Lift  to  drag  ratio 


C  =C 


D 


D 


est  /  esti 


c.  Density  gain 


Ka  =Ka 

Pest  Pestinil 


Figure  A.3:  Guidance  Initialization  Logical  Flow:  initial_guid.m 
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Table  A-4:  Inputs,  initial  guid.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

i 

IYD 

I  YD 

Unit  normal  to 
the  desired  orbit 
plane 

double[3] 

Derived 

i 

V_NAV 

V 

nav 

Vehicle  inertial 
velocity  vector 

ft/s 

double[3] 

Navigation 

Table  A-5:  Local  Guidance  Parameters,  initial_guid.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Use 

L 

COSPHIMIN 

C0S(^Lin 

Cosine  of  the 
maximum  bank 
angle 

■ 

double 

Initialized 

1 

COSPHIMAX 

C0S(^Lnx 

Cosine  of  the 
minimum  bank 
angle 

■ 

double 

Initialized 

1 

SIGNOFBANK 

Sign((f>) 

Sign  of  current 
bank  angle  guess 

U 

double 

Initialized 

1 

PHICMD 

AmcJ 

Commanded  bank 
angle 

deg 

double 

Initialized 

1 

SIGN_OF_ 

BANKJPAST 

Sisn(<f>)pasl 

Sign  of  previous 
bank  angle  guess 

■ 

double 

Initialized 

1 

COS_PHI_DES 

cos  (^J 

Cosine  of  the 
desired  bank 
angle 

■ 

double 

Initialized 

1 

CORRIDOR, 

SLOPE 

Aerocapture 
corridor  slope 

(-) 

double 

Initialized 

1 

LIFT_DOWN_ 

REVERSAL 

F 

1  LDR 

Flag  indicates 
need  for  a  lift 
down  reversal 

■ 

double 

Initialized 

1 

DPHI_ 

OVER_UNDER 

&OU 

Bank  angle  which 
gives  direction 
vehicle  takes  to 
reverse  direction 

deg 

double 

Initialized 
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for  a  sign  change 

1 

MODEL_ 

LIFTDOWN 

F 

1  MLD 

Flag  to 

activate/deactivate 
lift  down 
modeling 

1 

double 

Initialized 

1 

CPHI_ 

LIFT_DOWN 

Cosine  of  the  lift 
down  bank  angle 

double 

Initialized 

1 

SPHI_ 

LIFTDOWN 

sin  {<t>LD) 

Sine  of  the  lift 
down  bank  angle 

■ 

double 

Initialized 

l.AERO 

CDEST 

cD 

uest 

Drag  coefficient 
estimate 

■ 

double 

Initialized 

l.AERO 

L_OVER_D_EST 

L\ 

D)es, 

Lift  to  drag  ratio 
estimate 

double 

Initialized 

l.AERO 

KRHOEST 

KP„, 

Density  gain 
estimate 

double 

Initialized 

Table  A-6:  Constants,  initial_guid.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

c 

BANKMAX 

^max 

Maximum 
bank  angle 

deg 

double 

c 

BANK_MIN 

^min 

Minimum 
bank  angle 

deg 

double 

c 

DEG_TO_RAD 

c 

D-+R 

Degrees  to 
radians 
conversion 
factor 

rad/deg 

double 

c 

PHIEI 

■ 

Initial  bank 
angle  guess 
at  entry 
interface 

deg 

double 

c 

COSPHIDESINIT 

cos(4,L 

Initial 
cosine  of 
the  desired 
bank  angle 

(-) 

double 

c 

CORRIDORMAX 

^max 

Maximum 
lateral  out 

ft/s 

double 

Value 

165.0 

15.0 

n 

180.0 

90.0 

0.0 

2000.0 
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CORRIDOR  MIN 


CORRIDOR  V  MAX  I  r 


CORRIDOR  V  MIN 


DPHIOVER 
UNDER  PRED 


MODEL  LIFT  DOWN  I  F, 


PHI  LIFT  DOWN 


CDESTINITIAL  C 


L_OVER_D JNITIAL  L  ' 

Dy 


K  RHO  EST  INITIAL  K 


of  plane 
velocity 

Minimum 
lateral  out 
of  plane 
velocity 


Maximum 
forward 
velocity  for 
corridor 
slope 


Minimum 
forward 
velocity  for 
corridor 
slope 


Initial  bank 
angle 
which 
gives 

direction  to 
reverse 


Initial 
model  lift 
down  flag 


Lift  down 
bank  angle 


Initial  drag 
coefficient 
estimate 


Initial  L/D 
coefficient 
estimate 


Initial 

density 

gain 

coefficient 

estimate 


ft/s 

double 

ft/s 

double 

ft/s 

double 

deg 

double 

(-) 

double 

deg 

double 

(-) 

double 

(-) 

double 

(-) 

double 

115.0 


1.4286 
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A.3  Aerodynamic  Properties 

The  aero_properties.m  function  updates  the  aerodynamic  properties. 

The  logical  flow  of  the  Aerodynamic  Properties  algorithm  is  illustrated  in  Figure  A.4. 

The  data  inputs  to  aero_properties.m  are  described  in  Table  A-7. 

The  local  guidance  parameters  used  or  updated  in  aerojproperties.m  are  described  in 
Table  A- 8. 

The  constants  used  in  aero_properties.m  are  described  in  Table  A-9. 


1 .  Compute  various  unit  vectors 


vVR 


y bi  ~  Omh, 


ody 


ZDCL  ~ 


y_Bl  X  XVR 
|  y BI  X  XVR 


(A.5) 

(A.6) 

(A.7) 


Notice  that  assuming  zero  sideslip,  this  is  essentially  a  velocity  frame:  x  is  out  the 
‘nose’  in  the  velocity  direction  (it  differs  from  the  velocity  direction  by  the  angle  of 
attack),  y  is  out  the  ‘right  wing’,  and  z  completes  the  right  hand  coordinate  system 
and  is  generally  away  from  earth  (positive  up). 


2.  Compute  aerodynamic  accelerations.  Drag  is  in  the  negative  velocity  direction 
and  lift  is  perpendicular  to  the  velocity  direction  (positive  away  from  earth). 


_  —T  A 

adrag  ~  ~ansbXVR 


_  —T  « 

aiift  ~  ansbZDCL 


(A.  8) 
(A-9) 


3.  If  the  acceleration  due  to  drag  is  greater  than  zero,  execute  steps  a  through  g. 
Otherwise,  exit  the  procedure. 


a.  Convert  velocity  in  the  inertial  frame  to  velocity  in  the  body  frame 


vrel=QTmvr 


rel 


(A- 10) 


b.  Use  geometry  to  calculate  the  angle  of  attack 
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a  =  tan' 


kJA 

n„(il 


c.  Compute  the  estimated  drag  coefficient 

C,m=C,t+\f,m+C,  a}x 


(A.l  1) 


(A.  12) 


d.  Compute  the  standard  atmospheric  density  at  the  current  altitude  using  the 
Atmospheric  Model  procedure. 

e.  Estimate  the  atmospheric  density  at  the  current  altitude  from  the  measured 
drag 


Pmeas 


2 mEiadrag 


(A.  13) 


r  meas  \  —  \l 

CjJref  K\ 

f.  Compute  and  filter  the  density  bias  estimate  by  executing  steps  i  through  iii: 


Calculate  the  density  bias  estimate 


(A.  14) 


ii.  Ensure  the  density  bias  estimate  is  within  the  minimum  and  maximum 
values  by  taking  the  median  of  Kp ,  Kp^  ,  and  Kp^  . 

iii.  Filter  the  density  bias  to  get  the  new  estimate 

<a-15) 

g.  Compute  and  filter  the  lift  to  drag  ratio  estimate  by  executing  steps  i  through 


i.  Calculate  the  lift  to  drag  ratio  estimate 

L  _  ay, 

D  a,rna 


(A.  16) 


ii.  Ensure  the  density  bias  estimate  is  within  the  minimum  and  maximum 

values  by  taking  the  median  of  — ,  —  ,  and  —  I 

U  D  J  min  ^  J  max 

iii.  Filter  the  lift  to  drag  ratio  to  get  the  new  estimate 


K  A 


(A.  18) 
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Figure  A.4:  Aerodynamic  Properties  Logical  Flow:  aero_properties.m 
Table  A-7:  Inputs,  aerojproperties.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

i 

QJBTOI 

Qbi 

Vehicle  body- 
to-inertial 
quaternion 

(-) 

double[4] 

Navigation 

i 

A_NAV 

°nav 

Vehicle  inertial 
aerodynamic 
acceleration 
vector 

double  [3] 

Navigation 

i 

V_REL_MAG 

■ 

Vehicle  relative 
velocity  vector 
magnitude 

ft/s 

double 

Derived 

i 

VRELNAV 

v 

V  rel 

Vehicle  relative 
velocity  vector 

ft/s 

double  [3] 

Derived 

Table  A-8:  Local  Guidance  Parameters,  aero_properties.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Use 

l.AERO 

CDEST 

cD 

uest 

Drag  coefficient 
estimate 

■ 

Updated 

l.AERO 

LOVERDEST 

-1 

D)„ 

Lift  to  drag  ratio 
estimate 

V) 

double 

Updated 

l.AERO 

KRHOEST 

Kn 

Pest 

Density  gain 
estimate 

■ 

double 

Updated 

Table  A-9:  Constants,  aerojproperties.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

c 

YBODY 

yB 

Unit  vector 
out  ‘right 
wing’  of 
vehicle 

(-) 

double 

1 

c 

RADTODEG 

Cr->d 

Radians  to 
degrees 
conversion 

deg/rad 

double 

180.0 

n 

c 

CDZERO 

Parasite 

drag 

coefficient 

(-) 

double 

1.4286 

c 

CD  ALPHA 

Slope  of 
the  Cd  vs. 
a  graph 

(-) 

double 

0 

c 

CDALPHASQ 

Cd, 

a~ 

Slope  of 
the  Ca  vs. 
a2  graph 

(-) 

double 

0 

c 

MASSEI 

m  El 

Vehicle 
mass  at  El 

slugs 

double 

15.4783 

c 

S  REF 

$ ref 

■■ 

double 

12.163 

c 

K_RHO_MIN 

Kn 

Pm\n 

Min 

density 

bias 

(-) 

double 

0.5 
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c 

KRHOMAX 

Prrux 

Max 
density 
bias  est. 

(-) 

double 

2.0 

c 

K_RHOJFILTER_GAIN 

Ka 

Pfs 

Filter  gain 
for  density 
bias  est. 

(-) 

double 

0.05 

c 

LOVERDMIN 

B 

Min  L/D 
estimate 

(-) 

double 

0.56 

c 

LOVERDMAX 

L\ 

Dj 

max 

Max  L/D 
estimate 

(-) 

double 

0.14 

c 

L_OVER_D_FILTER_GAIN 

Filter  gain 
for  L/D 
estimate 

(-) 

double 

0.2 

A.4  Predictor-Corrector  Sequencer 

The  pcsequencer.m  function  alternately  runs  the  Predictor  and  Corrector  functions  to 
determine  the  desired  bank  angle. 


The  logical  flow  of  the  Predictor-Corrector  Sequencer  algorithm  is  illustrated  in  Figure 
A.5. 

The  data  inputs  to  pc_sequencer.m  are  described  in  Table  A- 10. 

The  local  guidance  parameters  used  or  updated  in  pc  sequencer.m  are  described  in  Table 
A-ll. 

The  constants  used  in  pc  sequencer.m  are  described  in  Table  A-12. 


1.  The  Predictor-Corrector  Sequencer  is  executed  at  a  specific  (user-defined) 
frequency  in  the  Aerocapture  Executive.  Each  time  it  is  executed,  it  begins 
counting  acceptable  tries  from  zero.  The  following  list  of  counting  and  storage 
variables  must  be  reinitialized  before  each  run  of  the  Predictor-Corrector 
Sequencer: 

Nhigh=0 
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N iow  —  0 

N,apl  =  0 
Ng00d  =0 


C0ScaP,  =°° 


COS 


(<U= 


00 


COS 


bracket 


=  00 


COS  ,  =  00 

ex  trap 


extrap 


:  00 


a  bracket  ~  00 


Aa-0 


Fb=0 


^)W 


=  0 


2.  If  velocity  magnitude  is  less  than  the  velocity  to  model  lift  down,  terminate  lift 
down  modeling  by  resetting  the  model  lift  down  flag  to  zero: 

Fmld  =  0 

3.  Select  the  apogee  criteria  by  determining  if  the  velocity  magnitude  is  greater  than 
the  velocity  for  the  tighter  apogee  criteria.  If  the  velocity  is  still  ‘fast’  the 
program  uses  a  looser  requirement  to  hit  the  target  apogee  because  at  faster 
velocities  the  profile  is  more  sensitive.  As  the  velocity  decreases,  the  apogee  miss 
must  be  less  in  order  to  hone  in  on  the  target. 

4.  Calculate  the  desired  bank  angle  by  executing  steps  a  through  g  a  user-defined 
number  of  times: 

a.  Calculate  the  bank  angle  to  try  and  the  cosine  of  the  bank  angle  to  try  using 
the  Corrector  procedure. 

b.  Calculated  the  predicted  apogee  using  the  Predictor  procedure. 

c.  Compute  the  apogee  miss 


&a  =  a  d  -a, 


t  arg  et 


(A.  19) 
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d.  If  the  apogee  miss  is  less  than  the  apogee  correct  criteria,  then  the  try  was 
acceptable.  Set  the  cosine  of  the  desired  bank  angle  equal  to  the  cosine  of  the 
bank  angle  tried  and  exit  the  procedure 

cos(^J  =  cos(^ 

e.  If  the  cosine  of  the  bank  angle  tried  is  greater  than  the  maximum  cosine  of 
bank  angle  allowed,  and  the  apogee  miss  is  positive,  (translated,  the  spacecraft 
is  asking  for  a  more  shallow  bank  angle  than  it  is  capable  of  and  still  missed 
high)  full  lift  down  is  required.  Set  the  cosine  of  the  desired  bank  angle  equal 
to  the  maximum  cosine  of  the  bank  angle  and  exit  the  procedure 

cos(^J=cos(^)max 

f.  If  the  cosine  of  the  bank  angle  tried  is  less  than  the  minimum  cosine  of  bank 
angle  allowed,  and  the  apogee  miss  is  negative,  (translated,  the  spacecraft  is 
asking  for  a  steeper  bank  than  it  is  capable  of  and  still  missed  low)  full  lift  up 
is  required.  Set  the  cosine  of  the  desired  bank  angle  equal  to  the  minimum 
cosine  of  the  bank  angle  and  exit  the  procedure 

cos(^J  =  cos(4™ 

g.  Determine  the  nature  of  the  solution  (i.e.  high,  low,  captured,  etc.)  and 
calculate  the  cosine  of  the  desired  bank  angle  by  executing  steps  i  through  vi: 

i.  If  the  prediction  captured  (determined  in  the  Predictor  procedure), 
increment  the  number  of  captured  solutions  by  one  and  set  the  capture 
cosine  value  equal  to  the  cosine  of  the  bank  angle  tried  and  proceed  to 
step  iii 

Ntv,  =  N,m+l  (A.20) 

cos„„,  =  cosfe, ) 

ii.  Otherwise,  save  the  Predictor  solution  as  a  ‘good’  solution  by 
executing  steps  1  through  4  and  then  proceed  to  step  iii: 

1 .  Increment  the  number  of  good  solutions 

<A-2I> 

2.  Rotate  the  new  good  solution  into  the  extrapolation  variables 

C0S extrap  (%)  =  C0S  extrap  (0 

cosei(ra/,(l)=cos(^) 


101 


®  extrap  (^  )  ^  extrap  ) 

® extrap  (0  ®  pred 

3.  If  the  predicted  apogee  is  higher  than  the  target  apogee,  the 
solution  is  high.  Increment  the  number  of  high  solutions  using 
Eq.  A. 22,  update  the  high  bracket  terms,  and  proceed  to  step  iii 

NUgh=Nhigh+\  (A.22) 

cosiracte,(l)  =  cos(^) 

® bracket  (0  ®  pred 

4.  Otherwise,  the  predicted  apogee  is  lower  than  the  target 
apogee,  and  the  solution  is  low.  Increment  the  number  of  low 
solutions  using  Eq.  A.23,  update  the  low  bracket  terms  and 
proceed  to  step  iii 

Nlow=Nlow  +  l  (A.23) 

c°sW£,(2)  =  cos(^) 

®  bracket  (^)  ^  precl 

iii.  If  there  is  a  high  and  a  low  (including  captured)  solution,  we  have 
bracketed  the  target  apogee.  Set  the  bracketed  flag  to  1 : 

Fb=  1 

otherwise,  proceed  to  step  iv. 

iv.  If  we  have  bracketed  the  target  apogee  execute  the  following  steps  and 
then  proceed  to  step  v 

1 .  Calculate  the  delta  bank  angle  from  the  last  try 

=  (A.24) 

2.  If  the  delta  bank  angle  is  less  than  the  minimum  allowable  delta 
bank  angle,  set  the  cosine  of  the  desired  bank  angle  to  the 
average  of  the  two  bank  angles 

cos(^J-cos  (A. 25) 

\  1  J 

otherwise,  proceed  directly  to  step  v. 

v.  Set  the  previous  bank  angle  tried  to  the  current  bank  angle 
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filry 

vi.  If  we  have  reached  the  maximum  number  of  runs  to  update  the  bank 
angle,  limit  the  predictions  by  executing  steps  1  and  2  and  then  exit 
the  procedure 

1 .  Correct  the  bank  angle  to  try  once  more  without  prediction  by 
running  the  Corrector  procedure. 

2.  Set  the  cosine  of  the  desired  bank  angle  equal  to  the  cosine  of 
this  most  recently  updated  bank  angle 

cos(^J  =  cos(^)min 
otherwise,  exit  the  procedure. 
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Figure  A.5:  Predictor-Corrector  Sequencer  Logical  Flow:  pc  sequencer.m 


Table  A-10:  Inputs,  pcsequencer.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

i 

V_NAV_MAG 

HI 

Vehicle  inertial 
velocity  vector 
magnitude 

ft/s 

double[3] 

Derived 

Table  A-ll:  Local  Guidance  Parameters,  pc_sequencer.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Use 

1 

COS_PHI_MIN 

cosMmi„ 

Cosine  of 
maximum 
bank  angle 

■ 

double 

Used 

1 

CO  SPHIM  AX 

C0S(^Lax 

Cosine  of 
minimum 
bank  angle 

■ 

double 

Used 

1 

COS_PHI_DES 

Cosine  of 
the  desired 
bank  angle 

■ 

double 

Updated 

1 

MODELLIFTDOWN 

F 

1  MLD 

Flag  which 
activates 
lift  down 
modeling 

1 

double 

Updated 
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Table  A-12:  Constants,  pc  sequencer.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

c 

DEGTORAD 

Degrees  to 
radians 
conversion 
factor 

Rad/deg 

double 

7C 

180.0 

c 

VIMODELLIFTDOWN 

vLD 

Velocity 
switch  for 
lift  down 
modeling 

ft/s 

double 

27900.0 

c 

VI_LOOSE_APOGEE 

vLA 

Switchover 
for  apogee 
tolerances 

ft/s 

double 

30000.0 

c 

APOGEEEPSILON 1 

ae\ 

Apogee 
tolerance 
limit  for 
fast 

velocities 

nm 

double 

25.0 

c 

APOGEEEPSILON2 

ae2 

Apogee 
tolerance 
limit  for 
slower 
velocities 

nm 

double 

1.0 

c 

APOGEETARGET 

^ target 

Target 

apogee 

nm 

double 

216.0 

c 

MAX_NUMBER_RUNS 

Maximum 
number  PC 
sequences 
(or  bank 
corrections) 
per  run 

(-) 

double 

4 

c 

INFINITY 

oo 

Infinity 

(-) 

double 

99999999 
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A.5  Corrector 


The  corrector.m  function  generates  a  new  bank  angle  to  try  in  Predictor. 

The  logical  flow  of  the  Corrector  algorithm  is  illustrated  in  Figure  A.6. 

The  data  inputs  to  corrector.m  are  described  in  Table  A-13. 

The  data  outputs  of  corrector.m  are  described  in  Table  A- 14. 

The  local  guidance  parameters  used  or  updated  in  corrector.m  are  described  in  Table 
A-15. 

The  constants  used  in  corrector.m  are  described  in  Table  A- 16. 

1 .  Determine  the  method  used  to  make  the  new  bank  angle  guess: 

a.  If  this  is  the  first  run,  use  Method  1 . 

b.  Otherwise,  if  we  have  bracketed  the  solution  and 

i.  The  number  of  low  solutions  is  not  zero,  we  have  a  high  solution 
and  a  low  solution,  use  Method  2. 

ii.  There  is  a  high  solution  and  a  capture,  or  more  than  one  high 
solution,  use  Method  3. 

c.  Otherwise,  if  we  have  not  bracketed  the  solution  and 

i.  There  is  only  one  good  solution,  use  Method  5. 

ii.  There  are  at  least  two  good  solutions,  use  Method  6. 

iii.  There  are  no  good  solutions  yet,  use  Method  4. 

2.  Execute  the  applicable  steps  for  each  method  to  determine  the  new  bank  angle 
guess: 

a.  Method  1 :  Run  the  last  guess  from  the  previous  guidance  cycle 

c°s(<U=cos(^J 

b.  Method  2:  Interpolate  a  high  guess  and  a  low  guess  to  the  target  apogee 

c.  Method  3:  Interpolate  a  high  guess  and  a  capture  or  interpolate  two  high 
guesses 
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i.  If  there  is  only  one  high  solution,  calculate  the  new  bank  angle 
guess  from  the  high  solution  and  the  captured  solution 

C0S(^rrv  )  ~  C0S  bracket  0)  +  {COScapt  ~  C0Sbracket  (0)Q.PC  (A. 26) 

ii.  Otherwise,  interpolate  a  new  bank  angle  to  try  from  the  two  good 
solutions.  If  this  new  bank  angle  to  try  is  greater  than  or  equal  to 
the  current  cosine  capture  value,  calculate  the  new  bank  angle 
guess  from  the  high  solution  and  the  captured  solution 

cos(^)=  cos^XO+lcos^-cos^XOferc  (A-27) 

d.  Method  4:  March  out  of  the  capture  region 

cos(^0, )  =  coscapl  -  (A.28) 

e.  Method  5:  Extrapolate  one  good  guess  using  a  stored  sensitivity 

i.  Calculate  sensitivity  if  the  velocity  magnitude  is  less  than  the 
current  break  value  (i) 

S  =  Css(ifnav\  +  Csi(i)  (A.29) 

ii.  Calculate  the  cosine  of  the  bank  angle  to  try 

cos  fe,, )  =  cos  \asi  +  y  (A.30) 

The  sensitivity-bank  angle  loop  can  run  up  to  6  times  (user 
defined).  The  final  bank  angle  guess  is  either  the  sixth  calculation 
or  the  (i- 1  )th  calculation  before  the  sensitivity  break  value  (i)  is  less 
than  the  velocity  magnitude. 

f.  Method  6:  Extrapolate  two  high  guesses  or  two  low  guesses  to  the  target 
apogee 

3.  Ensure  the  cosine  of  the  bank  angle  guess  is  within  the  minimum  and  maximum 

cosine  of  bank  angle  values  by  taking  the  median  of  cos(^(ry  )>  cosWmin’and 

cosM  • 

vr  'I  max 

4.  Calculate  the  bank  angle  to  try 

<t>try  =  C0S"'  C0Sfeo -  ))  (A’3 
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Figure  A.6:  Corrector  Logical  Flow:  corrector.m 
Table  A-13:  Inputs,  corrector.m 


inertial 

velocity 

magnitude 


Cosine  of 
bank  angle 
guess 


Flag  which 
indicates 
target 
bracketed 


COS  BRACKET 


Units 

Data 

Type/ 

Precision 

ft/s 

double[3] 

(-) 

double 

(-) 

double 

(-) 

double  [2] 

Source 
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bracket  the 
target 

APOGEEBRACKET 

^  bracket 

Apogees 
which 
bracket  the 
target 

nm 

double  [2] 

pc_sequencer 

CO  S_C  APT 

cos  wpr 

Cosine  of 
the  last  bank 
that  caused 
capture 

(-) 

double 

pcsequencer 

COS_EXTRAPOLATE 

cos„rru/, 

Cosine  of 
angles  for 
interpolation 

(-) 

double  [2] 

pc_sequencer 

APOGEEEXTRAPOLATE 

d  extrap 

Apogees  of 
solutions  for 
interpolation 

nm 

double  [2] 

pcsequencer 

NUMBER_GOOD 

^  good 

No.  good 
solutions 

(-) 

double 

pc  sequencer 

NUMBERHIGH 

N  high 

No.  high 
solutions 

(-) 

double 

pcsequencer 

NUMBERLOW 

Nlow 

No.  low 
solutions 

(-) 

double 

pc_sequencer 

DELTA  APOGEE 

A a 

Apogee 

error 

nm 

double 

pc_sequencer 

RUN  NUMBER 

nm 

Run  number 

(-) 

double 

pcsequencer 
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Table  A-14:  Outputs,  Corrector.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Destination 

PHI_TRY 

$  try 

Bank  angle 
guess 

deg 

double 

predictor 

COS_PHI_TRY 

Cosine  of  the 
bank  angle 
guess 

(-) 

double 

predictor 

Table  A-15:  Local  Guidance  Parameters,  corrector.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Use 

1 

COSPHIMIN 

cosWmi„ 

Cosine  of  the 
maximum 
allowable  bank 
angle 

double 

Used 

1 

COSPHIMAX 

C0S(^Lax 

Cosine  of  the 
minimum 
allowable  bank 
angle 

double 

Used 

1 

COSPHIDES 

C0S(^J 

Cosine  of  the 
desired  bank 
angle 

double 

Updated 
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Table  A-16:  Constants,  corrector.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

C 

RADTODEG 

r 

Radians  to 
degrees 
conversion 
factor 

deg/rad 

double 

180.0 

7C 

C 

LIFT_PERCENT_CAPTURE 

r 

^ LPC 

double 

0.50 

c 

LIFT_INC_CAPTURE 

Cue 

Cosine 
amount 
used  to 
march  out 
of  capture 
region 

double 

-0.10 

c 

SENBREAKPOINT 

Csbp 

Vector 
determines 
which 
stored 
sensitivity 
equations 
to  use 

double[6] 

34000.0 

32000.0 

30000.0 

27500.0 

26000.0 

0.0 

c 

SENJNTER 

Css 

y  intercept 
of  the 
stored 
sensitivity 
equation 

double[6] 

c 

SENSLOPE 

Csi 

Slope  of 
the  stored 
sensitivity 
equation 

double[5] 

'o  o  o  o  o  o' 

o  o  o  o  o  o 

1 _ 1 
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A.6  Predictor 


The  predictor.m  function  calculates  the  predicted  apogee  assuming  the  current  bank  angle 
guess  is  constant  for  the  remainder  of  the  trajectory.  It  also  determines  if  the  prediction 
captures. 

The  logical  flow  of  the  Predictor  algorithm  is  illustrated  in  Figure  A.7. 

The  data  inputs  to  predictor.m  are  described  in  Table  A- 17. 

The  data  outputs  of  predictor.m  are  described  in  Table  A- 18. 

The  local  guidance  parameters  used  or  updated  in  predictor.m  are  described  in  Table 
A- 19. 

The  constants  used  in  predictor.m  are  described  in  Table  A-20. 


1 .  Initialize  variables: 

a.  Predictor  state  vector 


b.  Calculate  predicted  altitude  (m  file  function) 

c.  Sine  and  cosine  of  the  predicted  bank  angle 

cos(^<J=cos(^J 

sin(^)=sin(^,) 

d.  Lift  to  drag  ratio 


(A.32) 

(A.33) 

(A.34) 

(A.35) 
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£ 

D 


J  pred 


£ 

DJest 


e.  Prediction  captured  flag  (a  value  of  1  indicates  the  prediction  captured  and  a 
value  of  zero  indicates  it  did  not) 


Fpc  —  0 

f.  Sign  of  the  sine  of  the  lift  down  bank  angle 


sin(^iZ) \.gn  =  sign(<f>)sm((/>LD  ) 

(A.36) 

Calculate  two  times  the  ballistic  coefficient 

„  „„  2  mF, 

2BC  = - — 

CD  S 

uesi 

(A.37) 

Execute  the  predictor  loop  for  a  user-defined  number  of  time  steps  by  performing 

steps  a  through  c: 

a.  Execute  the  Fourth  Order  Runge-Kutta  Integration  Loop  four  times  by 

performing  steps  i  through  vii: 

i.  Calculate  the  relative  velocity 

Kel  =  Kred  ~  l? pole  X  ^  pred  ] 

(A.38) 

1  y  I2  -  yT 

|  relpred  |  relpred 

(A.39) 

i_  j  /i—  |2 

|  F reiprcd  j  y  |  ^ 'relpred  \ 

(A.40) 

ii.  Use  the  Atmospheric  Model  procedure 

to  calculate  the  standard 

density  at  the  current  altitude. 

iii.  Compute  aerodynamic  accelerations 

F  p  P st d  [  F relpred  | 

°D  ~  BC 2 

(A-41) 

L' \ 

aL  ~  D  I  a° 

■D  s  pred 

(A.42) 

iv.  Compute  the  lift  vector 

v 

y  relpred 

vel  -  y 

relpred 

(A.43) 
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&* 

II 

X 

f 

S. 

(A.44) 

a 

II 

^1 

a 

^-Hl 

a 

(A.45) 

\hatc  |  | 

(A.46) 

7  ha, 

^  lot  It  1 

(A.47) 

latc 


1 .  If  lift  down  modeling  is  active,  and  the  velocity  is  less  than  the 
velocity  to  model  lift  down,  model  lift  down,  calculate  the  lift 
vector  using  the  model  lift  down  cosine  and  sine  values 

lift  =  {la,  X  he,  )cos {<f>,d  )+Ihlsign{sm{<f>ld  ))  (A.48) 


2.  Otherwise,  compute  direction  for  bank  angle  to  try 
l,/,  =  (la,  x  he,  )cos(^d  )+  4,  Sin(^ra, ) 

aaccel  -aihlif,  ~aohel 

v.  Compute  gravity  acceleration  with  J2 
-  R, 


TJ  = _ i 

^  lj? 


pred 


pred  | 

Z  - uT  P 

^  pred  u  predr9 


Upred  pred 


Ra 


R 


pred 


((l-5Zl,1'ppr,J+2Zpr,iP.) 


(A.49) 
(A.  50) 

(A.51) 
(A- 52) 

(A.53) 


a 


g 


—h^U. 

R  pred 


(A.54) 


vi. 


vii. 

viii. 


Compute  total  acceleration 

® pred  Aaccel 

Perform  Runge-Kutta  integration  using  the  Integrator  procedure 
Compute  state  parameters 

[  R  pred  Rpred 


(A.55) 


(A.56) 
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_ 

/i—  i2 

(A.57) 

|  R  pred  1  = 

y  \Rprcd  \ 

\v  I2  = 

pred  | 

VT  V 

pred  pred 

(A.58) 

|  ^ pred  | 

^ pred 

(A.59) 

b.  Check  for  atmospheric  exit  by  executing  the  following  steps: 

i.  Calculate  centrifugal  velocity  and  predicted  flight  path  angle 

VT  R 

D  _  pred  pred 

KPrcd  " 


R 


(A.60) 


pred 


f  R 


Y pred 


=  sin 


pred 


pred  \J 


(A-61) 


ii.  If  the  altitude  is  greater  than  the  altitude  at  which  we  are  still  in  the 
atmosphere,  and  the  centrifugal  velocity  is  positive,  stop 
integrating  and  proceed  to  step  4. 

c.  Check  for  atmospheric  capture  by  executing  the  following  steps: 


i.  Calculate  centrifugal  acceleration 


• 

(  \12 

^  pred  ^  pred 

r,,», 

^Vpred) \ 

^  pred 

(A.62) 


ii.  If  centrifugal  velocity  and  centrifugal  acceleration  are  both 
negative,  the  prediction  has  captured 

1 .  Reset  the  capture  flag  to  indicate  the  prediction  captured 

F  =1 

1  capt  1 

2.  Stop  integrating  and  proceed  to  step  4. 


4.  If  the  prediction  captured,  set  the  predicted  apogee  to  negative  infinity  and  exit 
the  procedure. 


a  pred  = 


5.  Otherwise,  compute  the  predicted  apogee  and  then  exit  the  procedure. 


P  pred 


A® 


(A.63) 
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Figure  A.7:  Predictor  Logical  Flow:  predictor.m 


Table  A-17:  Inputs,  predictor.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

i 

RNAV 

Rnav 

Vehicle  inertial 
position  vector 

ft 

double[3] 

Navigation 

i 

V_NAV 

Vnav 

Vehicle  inertial 
velocity  vector 

ft/s 

double[3] 

Navigation 

PHI  TRY 

$  try 

Bank  angle 
guess 

deg 

double 

Corrector 

Table  A-18:  Outputs,  predictor.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Destination 

APOGEEPREDICTED 

a  pred 

Predicted 

apogee 

nm 

double 

pc_sequencer 

PREDCAPTURE 

F 

capt 

Flag 
variable 
indicates  if 
prediction 
captured 

(-) 

double 

pc_sequencer 

Table  A-19:  Local  Guidance  Parameters,  predictor.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Use 

1 

SIGNOFBANK 

sign{<t>) 

Sign  of  the 
current  bank 
angle  guess 

■ 

double 

Used 

1 

MODELLIFTDOWN 

F 

1  MLD 

Flag  which 
activates  lift 
down 
modeling 

1 

double 

Used 

1 

CPHI_LIFT_DOWN 

Cosine  of  the 
lift  down 
bank  angle 

■ 

double 

Used 
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.AERO 


CD  EST 


Drag 

coefficient 

estimate 


double  Used 


double  I  Used 


Table  A-20:  Constants,  predictor.m 


Structure  Variable  Name  Symbol  Description  Units  Data 
Name  Type/ 


MASS  El 


S  REF 


DEG  TO  RAD  Cr 


EARTH  RATE  coa 


mass  in  El 
units 


Reference 

area 

Degrees  to 
radians 
conversion 


Angular 
rotation  rate 
of  the  Earth 


Unit  vector 
through  the 
Earth’s  pole 


Switch  for 
lift  down 
modeling 


Earth  J2 
constant 


Radius  of 
Earth 


Gravitational  ft3/sec2 
parameter  of 
Earth 


Entry 

interface 

altitude 


Data 

Type/ 

Precision 

Value 

double 

15.4783 

double 

12.163 

double 

n 

180.0 

double 

7.2921 1488e~5 

double[3] 

[-0.0042772340' 
-0.0000901672 
[  0.9999908485 

double 

27500.0 

double 

1082.63e-6 

double 

6378137 

0.3048 

double 

3.986005el4 

0.30483 

double 

400000.0 

c 

ALTCAPT 

^capt 

Altitude 
vehicle  is 
considered 
captured 

ft 

double 

200000.0 

c 

FT_TO_NM 

r 

^FT^NM 

Conversion 
factor  feet  to 
nautical 
miles 

Nm/ft 

double 

1.0 

6076.116 

c 

INFINITY 

00 

Infinity 

(-) 

double 

99999999 

A.7  Integrator 

The  integrator.m  function  uses  a  fourth  order  Runge-Kutta  method  to  integrate  the 
predicted  state  vectors  forward  one  time  step. 


The  logical  flow  of  the  Integrator  algorithm  is  illustrated  in  Figure  A.8. 
The  data  inputs  to  integrator.m  are  described  in  Table  A-21. 

The  data  outputs  from  integrator.m  are  described  in  Table  A-22. 

The  constants  used  in  integrator.m  are  described  in  Table  A-23. 


1 .  If  this  is  the  first  step  of  integration,  execute  the  following  steps: 

a.  Set  values  for  original  position  and  velocity 
=  R, 


orig 


'‘pred 


V  -v 

orig  pred 


b.  Set  accumulated  velocity  and  acceleration 


V  =  V 

accum  pred 


anrrum  =  a, 

accum  pred 


c.  Perform  first  step  of  integration 


pred 


pred  r  orig 


V  .  +-A ta 

f  nrifr  1  _  1 


pred 


2.  If  this  is  the  second  step  of  integration,  execute  the  following  steps: 


(A.  66) 
(A.67) 
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a.  Set  accumulated  velocity  and  acceleration 


-V  +2V 

accum  pred 

(A.68) 

&  accum  ^^pred 

(A.69) 

b.  Perform  second  step  of  integration  using  Eq  A. 66  and  Eq  A. 67 
3.  If  this  is  the  third  step  of  integration,  execute  the  following  steps: 

a.  Set  accumulated  velocity  and  acceleration  using  Eq.  A. 68  and  A.69 

b.  Perform  third  step  of  integration  using  Eq  A.66  and  Eq  A.67 


4.  Otherwise,  this  is  the  fourth  step,  perform  fourth  step  of  integration 


Figure  A.8:  Integrator  Logical  Flow:  integrator.m 
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Table  A-21:  Inputs,  integrator.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

INTEGLOOP 

F,l 

Count  variable  to 
determine  which 
step  of 
integration  to 
perform 

(-) 

double[3] 

predictor 

RPRED 

^ pred 

Predicted  vehicle 
inertial  position 
vector 

ft 

double 

predictor 

VPRED 

V 

pred 

Predicted  vehicle 
inertial  velocity 
vector 

ft/s 

double 

predictor 

APRED 

d  pred 

Predicted  vehicle 
acceleration 
vector 

"ft/s3"" 

double 

predictor 

Table  A-22:  Outputs,  integrator.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Destination 

RPRED 

Rprcd 

Predicted 
vehicle  inertial 
position  vector 

ft 

double 

predictor 

VPRED 

v 

pred 

Predicted 
vehicle  inertial 
velocity  vector 

ft/s 

double 

predictor 

Table  A-23:  Constants,  integrator.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

c 

DELT  ATPRED 

At 

Time  step  for 
integration 

s 

double 

20 
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A.8  Lateral  Guidance 


The  lateral  control.m  function  computes  the  sign  for  the  commanded  bank  angle  and 
determines  the  need  for  the  roll  reversals  that  keep  the  spacecraft  in  the  desired  orbital 
plane, 


The  logical  flow  of  the  Lateral  Guidance  algorithm  is  illustrated  in  Figure  A. 9. 

The  data  inputs  to  lateral_control.m  are  described  in  Table  A-24. 

The  local  guidance  parameters  used  or  updated  in  lateral_control.m  are  described  in 
Table  A-25. 

The  constants  used  in  lateral  control.m  are  described  in  Table  A-26. 


1.  Compute  the  maximum  lateral  corridor  out  of  plane  velocity  by  executing  the 
following  steps: 


a.  If  the  inertial  velocity  magnitude  is  less  than  or  equal  to  the  minimum  corridor 
velocity,  the  corridor  out  of  plane  velocity  is  equal  to  the  minimum  corridor 
out  of  plane  velocity 


C  =  C 

'“'min 

b.  Otherwise  calculate  the  maximum  corridor  out  of  plane  velocity  based  on  the 
current  forward  velocity  and  the  equation  of  the  line  and  ensure  it  is  within  the 
minimum  and  maximum  values 


C  =  Cmm  -sign 


£ 

D 


V... 


(A.72) 


est  J 


2.  Compute  the  actual  out  of  plane  velocity  based  on  the  navigational  velocity  and 
the  unit  normal  to  the  desired  plane 

y  =  Vjjyd  (A.73) 


3.  If  absolute  value  of  the  out  of  plane  velocity  is  greater  than  the  maximum  corridor 
out  of  plane  velocity  |j|  >  C,  the  vehicle  is  outside  of  the  corridor.  The  sign  of 
bank  is  opposite  that  of  the  out  of  plane  velocity  sign  to  cause  a  roll  reversal 

sign{<f>)  =  -sign(y)  (A.74) 
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4.  Otherwise,  the  sign  of  bank  is  equal  to  the  current  sign  of  bank  global  parameter 
value 

sign(<f>)  =  sign(<f)  (A.75) 


Figure  A.9:  Lateral  Guidance  Logical  Flow:  lateral_control.m 
Table  A-24:  Inputs,  lateral  control.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

i 

V_NAV_MAG 

■ 

Vehicle  inertial 
velocity  vector 
magnitude 

ft/s 

double 

Derived 

i 

V_NAV 

V 

nav 

Vehicle  inertial 
velocity  vector 

ft/s 

double[3] 

Navigation 

i 

IYD 

I  YD 

Unit  normal  to 
the  desired  orbit 
plane 

(-) 

double[3] 

Derived 
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Table  A-25:  Local  Guidance  Parameters,  lateralcontrol.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Use 

1 

CORRIDORSLOPE 

c 

^  corr 

Aerocapture 
corridor  slope 
(out  of  plane 
velocity  over 
forward 
velocity) 

(-) 

double 

Used 

1 

SIGNOFBANK 

sign(<t>) 

Sign  of  the 
current  bank 
angle  guess 

■ 

double 

Updated 

l.AERO 

LOVERDEST 

-1 

D)es< 

Lift  to  drag 
ratio  estimate 

■ 

double 

Used 

Table  A-26:  Constants,  lateral  control.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

C 

CORRIDORMAX 

^max 

Maximum 
lateral 
corridor  out 
of  plane 
velocity 

fit/s 

double 

2000.0 

C 

CORRIDOR_MIN 

^min 

Minimum 
lateral 
corridor  out 
of  plane 
velocity 

ft/s 

Double: 

15.0 

c 

CORRIDORVMAX 

Maximum 
forward 
velocity  to 
use  corridor 
slope 

ft/s 

double 

31000.0 
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A.9  Bank  Angle  Command  Incorporation 

The  cmdincorporation.m  function  updates  the  bank  angle  command. 

The  logical  flow  of  the  Bank  Angle  Command  Incorporation  algorithm  is  illustrated  in 
Figure  A.  10. 

The  data  inputs  to  cmd_incorporation.m  are  described  in  Table  A-27. 

The  data  outputs  from  cmd_incorporation.m  are  described  in  Table  A-28. 

The  local  guidance  parameters  used  or  updated  in  cmd_incorporation.m  are  described  in 
Table  A-29. 

The  constants  used  in  cmd  incorporation.m  are  described  in  Table  A-30. 

1 .  Select  the  nominal  bank  angle  command 

a.  If  inertial  velocity  is  less  than  the  bias  switch  velocity 

cos(^om)  =  cos(A/0i(l)) 

b.  Otherwise 

cos(^0J  =  cos(4OJ(2)) 

2.  Update  the  bank  command 

c.  If  inertial  velocity  is  less  than  the  command  bias  required  velocity,  bias  the 
commanded  bank  angle 

cos(^cmJ  =  cos(^es)+A:cA(cos(^)-cos(Aom))  (A.76) 

d.  Otherwise,  do  not  bias  the  commanded  bank  angle 

C°s(^mrf)=COS(^)  (A-77) 

3.  Limit  the  cosine  of  the  bank  angle 
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Figure  A.10:  Command  Incorporation  Logical  Flow:  cmd_incorporation.m 


Table  A-27:  Inputs,  cmd_incorporation.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

i 

V_NAV_MAG 

\V  1 

|  nav  | 

Vehicle  inertial 
velocity  vector 
magnitude 

ft/s 

double 

Derived 

Table  A-28:  Outputs,  cmdincorporation.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Destination 

PHICMD 

Commanded 
bank  angle 

deg 

double 

predguid 

Table  A-29:  Local  Guidance  Parameters,  cmd  incorporation.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Use 

1 

COSPHIMIN 

cosU]  . 

vr  'Imin 

Cos  max 
bank  angle 

■ 

double 

Used 

1 

COSPHIMAX 

C0S(^Lax 

Cos  min 
bank  angle 

double 

Used 
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Table  A-30:  Constants,  cnuHncorporation.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

c 

RADTODEG 

r 

^R->D 

Radians  to 
degrees 
conversion 
factor 

deg/rad 

double 

180.0 

n 

c 

VIBIASSWIT  CH 

Velocity  to 
switch 
between 
bias  angles 

ft/s 

double 

30800.0 

c 

COSPHIBIAS 

cos(4to 

Vector 
which 
stores  two 
bias  angles 

(-) 

double[2] 

g 

c 

VI_COMMAND_BIAS 

Minimum 
velocity  to 
bias  bank 
angle 

ft/s 

double 

26500.0 

c 

KCOMMANDBIAS 

Kcb 

Gain  which 
determines 
how 

quickly  the 
solution 
biases 

(-) 

double 

2 

A.10  Atmospheric  Model 
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The  atmos_model.m  function  computes  the  atmospheric  density  using  a  fourth  order 
curve  fit  of  scale  height. 

The  logical  flow  of  the  Atmospheric  Model  algorithm  is  illustrated  in  Figure  A.l  1. 

The  data  inputs  to  atmos_model.m  are  described  in  Table  A-3 1 . 

The  data  outputs  from  atmos_model.m  are  described  in  Table  A-32. 

The  constants  used  in  atmos  model.m  are  described  in  Table  A-3  3. 


1 .  Calculate  the  normalized  altitude 


= 


H 


H 


(A.78) 


ref 


2.  If  the  current  altitude  is  greater  than  the  maximum  atmospheric  altitude,  calculate 
the  scale  height  using  the  maximum  atmospheric  altitude 

HS„ 


TTO  __ 

1  ^  rt/irm 


alt  max 


H 


(A.79) 


ref 


3.  Otherwise,  use  a  fourth  order  curve  fit  to  get  scale  height 

HS„orm  =  C»sHlrm  +  C“sH\orm  +  C?Hlrm  +  Cf  Hnorm  +  Cx "s  (A.80) 

4.  Calculate  density 

Ps,d=Prefe  /mm  (A-81) 


FigureA.il:  Atmospheric  Model  Logical  Flow:  atmos_model.m 
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Table  A-31:  Inputs,  atmos_model.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Source 

H 

H 

Altitude 

ft 

double 

various 

Table  A-32:  Outputs,  atmos_model.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Destination 

RHOSTD 

Pstd 

Standard 

density 

double 

various 

Table  A-33:  Constants,  atmos_model.m 


Structure 

Name 

Variable  Name 

Symbol 

Description 

Units 

Data 

Type/ 

Precision 

Value 

c 

H_REF 

Href 

Reference 

altitude 

fit 

double 

242524.94 

c 

ATMOS_ALT_MAX 

^  alt  max 

Altitude  of 
maximum 
sensible 
atmosphere 

ft 

double 

400000.0 

c 

C_HS 

CHS 

Fourth 
order  curve 
fit  of 
density 

(-) 

double[2] 

“9.154583^-1" 

-2.574749e0 

3.023997e0 

-1. 59503  leO 
3.168398e-l_ 

c 

HS_ALT_MAX 

TT 

11  alt  max 

Maximum 

scale 

height 

ft 

i 

double 

20000.0 

c 

RHO  REF 

Pref 

Reference 

density 

9.9469294e-8 
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Appendix  B 
Simulation  Verification 

B.l  Dynamic  Simulation  Verification 

The  dynamic  simulation  written  in  Matlab  Simulink  was  verified  against  an  existing 
simulation.  The  existing  simulation  used  the  US  Standard  Atmosphere,  1976  and 
assumed  spherical  Earth.  These  conditions  were  applied  to  the  Simulink  simulation  and  a 
90°  constant  bank  angle  input  was  used.  Figure  B.l  illustrates  the  verification  setup 


Environm  ent 


Figure  B.l:  Dynamic  Simulation  Verification  Setup 

Plotting  both  the  trajectory  generated  by  the  above  set  up  and  the  trajectory  generated  by 
the  existing  simulation,  both  using  a  90°  constant  bank  angle  input,  Figure  B.2  through 
Figure  B.9  were  generated. 
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Figure  B.2:  Vehicle  Position  vs.  Time 
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Figure  B.3:  Vehicle  Velocity  Components  vs.  Time 


Figure  B.4:  Vehicle  Velocity  Magnitude  vs.  Time 


Figure  B.5:  Vehicle  Relative  Velocity  Components  vs.  Time 
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Figure  B.6:  Vehicle  Relative  Velocity  Magnitude  vs.  Time 
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Figure  B.7:  Vehicle  Acceleration  Components  vs.  Time 
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Figure  B.8:  Vehicle  Altitude  vs.  Time 
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Figure  B.9:  Vehicle  Acceleration  vs.  Time 


The  Simulink  simulation  was  then  updated  to  include  J2  effects  and  plotted  vs.  the 
verified  simulation  in  order  to  ensure  the  results  were  still  reasonable.  These  plots  are 
shown  in  Figure  B.10  through  Figure  B.17. 


x  107 


Figure  B.10:  Vehicle  Position  vs.  Time,  J2  Effects 
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Time 


Figure  B.ll:  Vehicle  Velocity  Components  vs.  Time,  J2  Effects 
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Figure  B.12:  Vehicle  Velocity  Magnitude  vs.  Time,  J2  Effects 


Figure  B.13:  Vehicle  Relative  Velocity  Components  vs.  Time,  J2  Effects 


Figure  B.14:  Vehicle  Relative  Velocity  Magnitude  vs.  Time,  J2  Effects 
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2.5 


Figure  B.17:  Vehicle  Acceleration  vs.  Time,  J2  Effects 

B.2  Guidance  Verification 

In  order  to  verify  the  guidance  portion  of  the  simulation  (PredGuid),  an  input  file  was 
used  to  run  the  Simulink  simulation  and  the  original  predguid  code  in  an  open  loop 
format.  This  setup  is  illustrated  in  Figure  B.18.  With  the  same  inputs  and  no  feedback, 
the  two  programs  generated  the  same  bank  angle  profile  as  shown  in  Figure  B.19. 


&-* 

Clock 


p  re  d  g  u  id 


Figure  B.18:  Open  Loop  Guidance  Verification  Setup 
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Figure  B.19:  Bank  Angle  Command  vs.  Time,  Open  Loop 

B.3  Closed  Loop  Simulation 

After  independently  verifying  the  dynamics  and  open  loop  guidance,  the  two  components 
were  run  in  a  closed  loop  format.  While  there  was  no  ‘known’  trajectory  to  verify 
against,  several  things  would  indicate  the  program  was  running  properly.  The  trajectory 
should  hover  around  90°  for  two  reasons:  First,  the  initial  conditions  used  in  the  closed 
loop  simulation  were  for  a  trajectory  that,  with  a  constant  bank  angle  of  90° ,  should  hit 
the  target  apogee.  Second,  the  command  biasing  angle  was  set  to  90° .  As  a  result,  the 
commanded  bank  angle  should  settle  at  90°.  The  closed  loop  setup  is  illustrated  in 
Figure  B.20.  The  resulting  commanded  bank  angle  profile  is  shown  in  Figure  B.21 . 


Memory  Block 


Figure  B.20:  Closed  Loop  Guidance  Verification  Setup 
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Time  (s) 

Figure  B.21:  Bank  Angle  Command  vs.  Time,  Closed  Loop 

This  profile  was  not  as  close  to  settling  at  90°  as  hoped.  The  problem  was  determined  to 
be  that  guidance  used  a  different  atmosphere  model  than  the  dynamics  were  using.  In 
reality,  guidance  does  not  have  perfect  knowledge  of  the  atmosphere  and  the  values  used 
by  guidance  would  not  be  identical  to  the  actual  environment.  However,  for  debugging 
and  development  purposes,  assuming  perfect  knowledge  of  the  atmosphere  helps  to 
determine  if  the  guidance  logic  itself  is  flawed  versus  not  being  robust  enough  to  handle 
atmospheric  dispersions.  As  a  result,  the  atmosphere  model  in  the  guidance  was  replaced 
with  the  US  Standard  Atmosphere,  1976.  Figure  B.22  shows  that  the  results  were  much 
closer  to  what  was  expected. 
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Figure  B.22:  Bank  Angle  Command  vs.  Time,  Closed  Loop,  US  1976 


B.4  Guidance  Corrections 


In  the  process  of  translating  the  original  guidance  code  from  Matlab  to  Simulink,  several 
errors  and  problems  were  discovered  in  the  code.  First,  the  simulation  was  not  executing 
bank  reversals.  This  lack  of  reversals  is  obvious  from  Figure  B.22.  However,  a  plot  of 
the  lateral  corridor  and  out-of-plane  velocity  (Figure  B.23)  further  verified  this 
observation.  As  the  vehicle  approached  the  lateral  corridor  limits,  instead  of  reversing 
bank  to  stay  within  the  limits,  it  continued  to  move  out  of  plane. 

Examination  of  the  code  revealed  that  the  unit  normal  to  the  desired  plane  was  being 
calculated  incorrectly.  Eq.  B.l  shows  the  correct  calculation  and  B.2  shows  the  original 
incorrect  calculation. 


vxr 

vxr 


(B.l) 


r  x  v 
!r  xv 


(B.2) 
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Figure  B.24  shows  the  lateral  corridor  and  out  of  plane  velocity  after  making  this 
correction. 


Figure  B.23:  Lateral  Corridor  and  Out  of  Plane  Velocity,  No  Bank  Reversals 


Figure  B.24:  Lateral  Corridor  and  Out  of  Plane  Velocity,  Corrected 


Further  examination  of  the  code  revealed  several  problems  with  sign  usage.  When  the 
code  was  originally  written  for  the  Aeroassist  Flight  Experiment  (AFE),  a  negative  value 
for  lift-to-drag  ratio  was  used.  When  translated  to  Matlab  for  another  aerocapture 
proposal,  the  decision  was  made  to  use  a  more  intuitive,  positive  lift-to-drag  ratio.  In  the 
process  of  translating  and  switching  signs,  several  instances  of  incorrect  sign  usage  were 
overlooked.  After  correcting  all  of  these  instances,  a  closed  loop  simulation  was  run  and 
produced  the  bank  angle  profile  in  Figure  B.25. 


Figure  B.25:  Bank  Angle  Command  vs.  Time,  Correct  Guidance 
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